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ABSTRACT 


"  Formulation  and  application  of  a  windward  plane  of  symmetry  laminar,  transitional, 
and  turbulent  boundary-layer  analysis  is  presented  for  sharp  cones  at  incidence  in  a 
supersonic  or  hypersonic  flow.  The  governing  boundary-layer  equations  in  the  plane  of 
symmetry  are  numerically  integrated  on  a  digital  computer  using  a  marching  implicit 
finite-difference  technique.  A  so-called  invariant  model  of  turbulence  is  used  in  a  two-layer 
eddy  viscosity-mixing  length  approach  for  calculation  of  the  turbulent  boundary  layer  in 
conjunction  with  an  intermittency  factor  treatment  of  the  transition  zone.  Comparison 
of  the  present  theory  with  experimental  data  (surface  heat  transfer,  boundary-layer 
parameters  such  as  displacement  thickness,  and  boundary-layer  profiles)  under  supersonic 
and  hypersonic  flow  conditions  over  sharp  cones  at  incidence  reveals  good  agreement. 
In  general,  smaller  crossflow  (outflow)  effects  on  the  windward  plane  of  symmetry 
boundary  layer  can  be  expected  for  turbulent  layers  as  compared  with  laminar  layers  subject 
to  the  same  boundary  conditions. 
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SECTION  I 
INTRODUCTION 

The  laminar,  transitional,  and  turbulent  boundary  layer  along  the  windward  ray  of 
a  sharp  cone  at  incidence  is  of  practical  importance  in  many  aeronautical  applications. 
An  excellent  example  is  the  compression  inlet  for  high-performance  engines  on  supersonic 
and  hypersonic  air-breathing  vehicles  at  angle  of  attack.  Three-dimensional  effects  on 
boundary-layer  behavior  are  also  of  interest  relative  to  maneuverable  reentry  bodies  and 
lifting  surfaces  (such  as  delta  wings)  at  incidence.  All  of  these  flows  involve  pressure 
gradients  in  the  plane  of  symmetry,  either  favorable  or  adverse,  and  crossflow  pressure 
gradients  that  may  cause  either  inflow  toward  or  outflow  from  the  plane  of  symmetry. 
An  excellent  collection  of  current  topics  in  three-dimensional  inviscid  and  viscous  analysis 
techniques  for  such  flows  as  discussed  above  may  be  found  in  Ref.  1. 

The  present  report  concentrates  on  formulation  of  a  windward  plane  of  symmetry 
laminar,  transitional,  and  turbulent  boundary-layer  analysis  for  sharp  cones  at  incidence 
in  a  supersonic  or  hypersonic  flow  and  comparison  of  results  with  those  from  experiments. 
The  governing  boundary-layer  equations  in  the  plane  of  symmetry  have  been  numerically 
integrated  on  a  digital  computer  using  a  marching  implicit  finite-difference  technique,  with 
a  two-layer  eddy  viscosity-mixing  length  approach  based  on  a  so-called  invariant  model 
of  three-dimensional  turbulence  used  for  calculation  of  the  turbulent  boundary  layer.  The 
transition  region  is  treated  through  an  intermittency  factor  approach  coupled  to  the  eddy 
viscosity  model  mentioned  above.  In  order  to  help  understand  the  basis  for  the  current 
analysis,  a  short  review  of  applicable  boundary-layer  technology  relative  to  analysis  of 
fully  turbulent  two-dimensional  flows  (utilizing  accurate  numerical  integration  techniques 
in  conjunction  with  eddy  viscosity  models  of  turbulence)  as  well  as  laminar  and  turbulent 
plane  of  symmetry  applications  is  given  below.  Also  included  is  a  summary  of  very  recent 
analytical  work  involving  numerical  calculation  of  transitional  boundary  layers. 

1.1  REVIEW  OF  BOUNDARY-LAYER  TECHNOLOGY  APPLICABLE  TO  LAMINAR, 

TRANSITIONAL,  AND  TURBULENT  PLANE  OF  SYMMETRY  FLOWS 

Calculation  techniques  for  the  two-dimensional,  nonsimilar,  fully  turbulent,  supersonic 
and  hypersonic  boundary  layer  utilizing  large  digital  computer  codes  have  become  available 
during  the  past  few  years,  e.g.,  Smith  and  Cebeci  (Ref.  2),  Patankar  and  Spalding  (Ref. 
3),  Bushnell  and  Beckwith  (Ref.  4),  Herring  and  Mellor  (Ref.  5),  Martellucci,  Rie,  and 
Sontowski  (Ref.  6),  Mayne  and  Dyer  (Ref.  7),  Cebeci  (Ref.  8),  and  Pletcher  (Ref.  9). 
Comparison  of  results  from  the  above  analyses  with  experimental  data  has  shown  good 
agreement  for  fully  turbulent  supersonic  and  hypersonic  (up  to  M„  ~  1 0)  boundary  layers 
on  sharp  and  slightly  blunted  cones,  as  well  as  sharp  flat  plates.  These  results  suggest 
that  the  eddy  viscosity  model  employing  the  Prandtl  mixing  length  is  indeed  adequate 
for  prediction  of  hypersonic  turbulent  boundary-layer  flows  over  simple  bodies  under  cold 
wall  conditions  (and  may  be  used  with  some  confidence). 

With  regard  to  calculations  for  transitional  boundary-layer  flows,  little  has  been 
accomplished.  Most  of  the  above-listed  references  employ  an  instantaneous  step  transition 
from  fully  laminar  to  fully  turbulent  flow  at  either  a  specified  body  location  or  local 
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Reynolds  number  along  the  boundary-layer  edge.  In  order  to  provide  a  suitable  progressive 
transition  model  for  the  transition  region,  Harris  (Ref.  10)  and  Adams  (Ref.  11)  have 
applied  an  intermittency-factor,  phenomenological  description  of  the  transition  region  based 
on  the  experimentally  observed  fact  that  transition  in  a  boundary  layer  is  characterized 
by  the  intermittent  appearance  of  turbulent  spots  which  grow  and  move  downstream  with 
the  fluid.  The  basis  behind  this  approach  may  be  found  in  the  work  of  Masaki  and  Yakura 
(Ref.  12).  In  general,  the  analyses  by  Harris  and  Adams  indicate  that  the  intermittency 
factor  treatment  of  the  transition  region  yields  a  reasonably  accurate  description  of 
transitional  heating  on  sharp  plates  and  cones  under  supersonic  and  hypersonic  conditions. 

The  mathematical  theory  of  the  three-dimensional  laminar  boundary  layer  as 
formulated  by  Moore  (Ref.  13)  and  Hayes  (Ref.  14)  has  been  available  for  about  twenty 
years.  Both  Moore  and  Hayes  showed  that  the  laminar  boundary-layer  equations  for  a 
sharp  cone  at  incidence  allow  a  similarity  transformation  which  reduces  the  number  of 
independent  variables  from  three  to  two.  Moore  (Ref.  15)  published  a  solution  to  these 
equations  in  the  limit  of  very  small  angle  of  attack,  corresponding  to  a  perturbation  solution 
for  the  invisdd  flow  at  small  angle  of  attack.  Moore  expressed  the  boundary-layer  profiles 
as  the  sum  of  a  basic  flow  (satisfying  the  zero  angle-of-attack  equations)  and  correction 
terms  proportional  to  angle  of  attack.  He  reduced  the  correction  terms  to  universal 
functions;  however,  his  results  were  valid  only  for  the  adiabatic  wall  case.  In  1952,  Moore 
(Ref.  16)  published  a  study  of  the  conical  boundary-layer  equations  at  large  angle  of 
attack.  In  this  work,  attention  was  restricted  to  the  symmetry  planes  where  the  equations 
reduce  to  a  set  of  ordinary  differential  equations.  Again,  he  considered  only  the  case 
of  zero  heat  transfer  to  the  surface,  although  extensions  to  include  heat  transfer  along 
the  windward  plane  of  symmetry  were  later  published  by  Reshotko  (Ref.  1 7)  and  Brunk 
(Ref.  18).  In  Ref.  16,  Moore  obtained  numerical  solutions  to  the  equations  on  the  windward 
symmetry  plane  for  all  angles  of  attack  and  on  the  leeward  symmetry  plane  for  very 
small  angles  of  attack.  A  recent  report  by  Bodonyi  and  Reshotko  (Ref.  19)  generalizes 
the  small  angle-of-attack  solutions  obtained  by  Moore  (Ref.  1 5)  by  removing  the  restrictions 
of  insulated  surfaces  and  unit  Prandtl  number. 

Recent  studies  by  Libby,  et  al.  (Ref.  20)  and  Trella  and  Libby  (Ref.  21)  have  been 
concerned  with  similar  solutions  for  the  hypersonic  laminar  boundary  layer  near  a  plane 
of  symmetry  allowing  surface  heat  transfer.  These  analyses  of  similar  flows  lead  to  the 
usual  restrictions  that  arise  for  similar  flows  without  crossflow  and  to  additional  restrictions 
associated  with  the  three-dimensionality.  A  linearization  procedure  in  conjunction  with 
an  integral  method  is  followed  to  describe  weak  three-dimensional  effects  in  Ref.  20;  a 
numerical  solution  of  the  full  nonlinear  system  allowing  arbitrary  favorable  and  adverse 
pressure  gradients  and  for  inflow  and  outflow  is  presented  in  Ref.  21.  The  manner  of 
numerical  solution  has  been  further  refined  for  the  plane  of  symmetry  problem  by  Libby 
and  Liu  (Ref.  22)  based  on  the  use  of  quasilinearization.  A  recent  paper  by  Libby  (Ref. 
23)  examines  the  three-dimensional  laminar  boundary  layer  near  the  windward  generator 
of  a  sharp  cone  at  angle  of  attack  subjected  to  uniform  mass  transfer,  either  suction  or 
injection. 

The  problem  of  turbulent  boundary-layer  flow  along  a  plane  of  symmetry  has  received 
little  attention  in  the  literature.  The  general  three-dimensional  turbulent  boundary-layer 


2 


AEDC-TR-71  -235 


equations  have  been  derived  by  Vaglio-Laurin  (Ref.  24)  and  Braun  (Ref.  25).  For  the 
case  of  incompressible  plane  of  symmetry  turbulent  boundary  layers,  Johnston  (Ref.  26) 
applied  an  integral  technique  yielding  results  which  agreed  well  with  experiment.  Mellor 
(Ref.  27)  reconsidered  the  same  incompressible  case  using  an  exact  numerical  integration 
of  the  governing  plane  of  symmetry  equations  in  conjunction  with  an  isotropic  model 
of  three-dimensional  turbulence  through  an  eddy  viscosity  approach.  Again,  good  agreement 
with  the  data  used  by  Johnston  was  obtained  which  supported  the  assumption  that  the 
eddy  viscosity  was  isotropic.  Very  recently  Zakkay  and  Calarese  (Ref.  28)  have  examined 
the  compressible  turbulent  boundary  layer  along  a  plane  of  symmetry  undergoing  both 
adverse  pressure  gradient  and  crossflow  effects.  The  three-dimensional  compressible  integral 
equations  are  written  along  the  symmetry  plane  and  integrated  numerically  using  a 
power-law  assumption  for  the  streamwise  velocity  profile  and  a  modified  Crocco  relation 
for  the  enthalpy  profile.  Reasonable  results  relative  to  experimental  data  for  boundary-layer 
thickness,  displacement  thickness,  and  momentum  thickness  were  reported  for  the 
windward  side  of  the  meridian  plane  along  a  sharp  cone  in  a  supersonic  stream. 

To  the  present  author's  knowledge,  the  current  investigation  is  the  first  to  consider 
the  plane  of  symmetry  problem  allowing  for  laminar,  transitional,  and  turbulent  flow  in 
a  single  unified  analysis  based  on  accurate  numerical  integration  of  the  governing 
boundary-layer  equations. 


SECTION  II 

ANALYTICAL  ANALYSIS 

The  present  analytical  investigation  employs  a  three-dimensional  windward  plane  of 
symmetry  laminar,  transitional,  and  turbulent  boundary-layer  analysis  coupled  with  a 
three-dimensional  invisdd  conical  flow  analysis  for  a  sharp  cone  at  incidence  in  a  supersonic 
or  hypersonic  stream.  Full  development  of  the  boundary-layer  analysis  will  be  presented 
below;  the  invisdd  analysis  utilizes  a  documented  digital  computer  code  which  will  be 
described  briefly. 

2.1  GOVERNING  BOUNDARY-LAYER  EQUATIONS 

The  present  analysis  employs  the  three-dimensional  compressible  turbulent 
boundary-layer  equations  in  terms  of  time-averaged  mean  flow  quantities  as  derived  by 
Vaglio-Laurin  (Ref.  24).  The  coordinate  system  is  taken  to  consist  of  geodesics  and  geodesic 
parallels  following  Moore  (Ref.  13).  In  this  coordinate  system  the  body  surface  is  defined 
by  y  =  0  and  a  point  in  space  is  defined  by  the  distances,  x,  y,  and  r(x)^  where  the 
length  r(x)$  depends  implicitly  on  the  distance  x  and  where  r(x)  has  the  dimensions  of 
length.  This  coordinate  system  is  especially  useful  in  the  analysis  of  flow  about  bodies 
for  which  a  coordinate  x  can  be  defined  such  that  body  cross  sections  are  similar  for 
various  values  of  x.  The  quantity  r(x)  then  gives  the  variation  of  scale  for  these  cross 
sections.  For  the  sharp  cone  geometry  under  present  consideration,  the  geodesic  coordinates 
are  taken  to  be  the  cone  generators  and  the  geodesic  parallels  are  the  circles  swept  by 
the  meridianal  angle.  The  corresponding  length  function  r(x)  is  the  local  radius  of  the 
body.  See  Fig.  1  for  clarification  of  the  sharp  cone  geometry,  nomenclature,  and  coordinate 
system.  The  velocity  components  are  taken  to  be  u,  v,  and  w  in  the  directions  of  x, 
y,  and  0,  respectively.  The  governing  equations  of  motion  are,  following  Vaglio-Laurin 
(Ref.  24): 
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CONTINUITY 
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where 


V  =  v+  4^-  (6) 

P 

and  the  usual  expressions  for  the  mean  and  fluctuating  parts  of  the  dependent  variables 
are  used;  e.g., 


P  =  P  +  P'  (7) 

Implicit  in  Vaglio-Laurin's  derivation  of  the  above  equations  are  the  following  stipulations: 

a.  The  rates  of  change  of  the  mean  flow  properties  in  the  x  and  <f>  directions 
[0(1>]  are  smaller  than  the  rates  of  change  in  the  y  direction  IC(8'1)]  by 
an  order  of  magnitude. 

b.  Mean  squares  and  products  of  the  turbulent  fluctuations  are'0(5);  that  is, 
the  turbulent  level  is  small.  The  terms  involving  mean  squares  of  the  velocity 
fluctuations  are  taken  to  be  negligible,  which  is  valid  for  high  Reynolds 
number  flows  with  a  zero  or  favorable  pressure  gradient. 

c.  The  time  average  molecular  transports  are  approximated  by  those  pertaining 
to  the  mean  flow  properties;  indeed,  even  the  latter  are  negligible,  except 
very  near  the  wall,  compared  with  terms  involving  the  turbulent  transports. 
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The  energy  equation  (5)  is  defined  in  terms  of  the  mean  stagnation  enthalpy 


H  =  h 


— 2  —  2 

UZ+  W* 


(8) 


If  the  two  momentum  equations  (2)  and  (3)  are  multiplied  by  u  and  w,  respectively, 
and  added,  one  obtains  an  equation  for  (u2  +  w2)/2  which  is  simplified  by  the  fact  that 
the  curvature  terms  (those  involving  3r/3x)  vanish.  If  then  this  equation  is  combined  with 
the  energy  equation  (5),  one  obtains  the  following  energy  equation  in  terms  of  the  mean 
static  enthalpy: 

ENERGY 


__5h 
pu  37 


—  ydh 


p  w  <5h 

r  dd> 


—  dp  w  dp 
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\37/  +W. 
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57  Fr37  "  p 
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(9) 


This  is  the  form  of  the  energy  equation  used  in  the  present  analysis. 

If  subscript  w  denotes  wall  and  subscript  e  denotes  outer  edge  of  the  boundary  layer, 
the  associated  boundary  conditions  on  the  above  defined  equations  are 

MOMENTUM 


ENERGY 


:u  =  v=  w=  uv  =  v  w  -  p  \  =  0 

lim  y-»~  :  u"  ->  U  e,  w  ->  We 

u  V '  -*  0,  v'w'-*b,  p  v '  -*  0 

y  =  0:h  =  hw,7T7=0 

lim  y  -*  oo  :  h  -»  hg,  v'h  ' 0 


(10) 


(11) 


which  reflect  the  requirements  of  no  slip  or  homogeneous  mass  injection  (suction  or 
blowing)  at  the  wall  as  well  as  a  prescribed  wall  enthalpy.  The  normal  momentum  equation 
(4)  reveals  that  the  static-pressure  variation  across  the  boundary  layer  is  negligible,  and 
hence  the  static  pressure,  p (x,0),  is  regarded  as  external  input  to  the  boundary-layer  analysis 
from  a  separate  inviscid  analysis.  The  outer  edge  velocities,  Ue  and  We,  as  well  as  the 
outer  edge  static  enthalpy,  he,  must  be  determined  from  the  inviscid  analysis  consistent 
with  the  imposed  static-pressure  distribution. 


5 


AEDC-TR-71-235 


The  gas  model  adopted  for  the  present  study  is  thermally  and  calorically  perfect 
air  having  a  constant  specific  heat  ratio  7  =  1.40  and  obeying  the  equation  of  state 

P  =  pRT  02) 

where  R  =  1716  ft2/sec2-°R.  Hence,  under  this  assumption, 


h  =  CpT  03) 

where  Cp  =  6006  ft2/sec2-°R.  The  laminar  viscosity,  n,  is  taken  to  obey  Sutherland's 
law 


JL- 


^ref 


Tref  +  198.6 
T  +  198.6 


(14) 


where  T  must  have  units  of  °R  and  subscript  ref  denotes  a  reference  condition.  The  laminar 
Prandtl  number,  Pr,  is  taken  to  be  a  constant  value  of  0.71  across  the  entire  boundary 
layer. 


2.2  TURBULENT  TRANSPORT  LAWS 


Before  Eqs.  (1),  (2),  (3),  and  (9)  can  be  solved,  expressions  must  be  supplied  for 
the  Reynolds  stress  or  turbulent  shear  terms  in  the  momentum  equations  and  the  turbulent 
flux  of  static  enthalpy  in  the  energy  equation.  The  approach  used  in  the  present  analysis 
is  to  model  these  terms  as  functions  of  the  mean-flow  variables  following  the  analysis 
by  Hunt,  Bushnell,  and  Beckwith  (Ref.  29). 

The  concept  that  the  Reynolds  stress  in  turbulent  flow  is  proportional  to  a  momentum 
exchange  coefficient  times  the  mean-flow  velocity  gradient  normal  to  the  surface  is  well 
known  and  commonly  used  in  turbulent  boundary-layer  analyses.  This  concept  is  based 
on  an  assumed  analogy  between  the  so-called  eddy  viscosity  and  the  molecular  viscosity. 
The  total  shear  components  in  the  streamwise  (x)  and  circumferential  (0)  directions  are 
written  as 


T 


X 


™  /  / 
-/JUV 


du  du 

+  '*dj 


(15) 


V  = 


p  v'w' 


<?w  <?w 

*W+t+W 


(16) 


where  the  eddy  viscosities  ex  and  in  the  x-  and  0-directions,  respectively,  might  in 
general  be  different.  Since  the  total  resultant  shear  must  be  a  vector  quantity,  its  magnitude 
is  written  as 


1/2 


(17) 
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In  a  similar  manner  the  total  heat  flux  in  the  energy  equation  (9)  becomes 


q  = 


db 

<?v  P 


■'h' 


k  <?h  k  db 

C  <Jv  +  C  <9v 

p  p  J 


(18) 


where  tc  is  the  so-called  eddy  thermal  conductivity. 

2.3  EDDY  VISCOSITY  MODEL 


The  simplest  approach  to  the  formulation  of  models  for  the  Reynolds  stress  is  based 
on  Prandtfs  mixing-length  hypothesis.  For  two-dimensional  flow,  this  hypothesis  states 
that 


U  » 

s 

V  « 


The  turbulent  shear  and  eddy  viscosity  are  then 


from  whence 


du 


fturb  -  "Pu'v'=  %=  Pi 
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du 
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(20) 


«  = 


(21) 


The  quantity  2,  is  termed  the  mixing-length  and  is  some  characteristic  length  related  to 
the  size  or  scales  of  eddies  responsible  for  the  flux  of  momentum  in  the  y-direction. 
Although  the  details  of  such  a  transfer  mechanism  are  not  well  understood,  the  basic 
concept  gives  satisfactory  results  even  at  hypersonic  Mach  numbers  in  the  presence  of 
heat-transfer  and  pressure  gradients  as  mentioned  in  the  review  of  turbulent  boundary-layer 
literature  given  in  Section  I. 


The  turbulent  boundary-layer  studies  by  Prandtl  (Ref.  30)  and  Glushko  (Ref.  31) 
are  based  on  the  concept  that  the  eddy  viscosity  should  depend  only  on  the  properties 
of  the  turbulence  and  a  local  eddy  scale.  Extending  this  concept  to  the  three-dimensional 
case  of  present  interest  suggests  that  the  eddy  viscosity  should  be  a  scalar  function 
independent  of  the  coordinate  direction.  Accordingly,  the  components  of  the  Reynolds 
stress  are  written  as 


'turb. 
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so  that  the  eddy  viscosity  is 


€ 


=  € 


O 


— p2  <?G 

■ 


(24) 


To  determine  the  scalar  function  G,  Eq.  (24)  is  substituted  into  Eq.  (17)  with  the  result 
that  the  magnitude  of  the  total  resultant  shear  becomes 


(25) 


Now,  by  analogy  with  Eqs.  (20)  and  (21)  where  the  velocity-gradient  function  for  the 
turbulent  shear  component  is  repeated  in  the  eddy  viscosity  expression,  the  relationship 
for  the  determination  of  the  G  function  becomes 


dG 

dy 


+ 


(26) 


which  results  in  the  final  form  of  the  eddy  viscosity 


(27) 


Hunt,  Bushnell,  and  Beckwith  (Ref.  29)  call  the  above  representation  of  the  Reynolds 
stress  the  invariant  turbulence  model.  The  key  assumption  in  this  model  is  that  the  eddy 
viscosity  is  a  scalar  function  independent  of  coordinate  direction.  In  a  recent  analysis 
of  three-dimensional  incompressible  turbulent  boundary-layer  flows  using  the 
kinetic-energy-of- turbulence  approach,  Nash  (Ref.  32)  has  advanced  arguments  that  the 
turbulent  shear  stress  is  likely  to  act  in  the  mean  rate  of  strain  direction,  defined  by 
the  components  of  the  mean  velocity  gradient  vector,  so  that  his  closure  equation  is  simply, 
in  the  nomenclature  of  the  present  analysis, 

rturb,  x  rturb,  tf> 

5u  ~  (28) 

<3y  dy 

The  same  equation  results  from  the  scalar  eddy  viscosity  model  presented  by  Eqs.  (22) 
and  (23)  above.  On  the  other  hand,  Bradshaw  (Ref.  33)  has  derived  a  set  of  differential 
equations  for  the  two  components  of  the  turbulent  shear  stress  based  again  on  the 
kinetic-energy-of-turbulence  approach  but  permitting  the  turbulent  shear  stress  vector  to 
deviate  from  the  mean  rate  of  strain  direction.  However,  in  the  near-wall  region,  Bradshaw's 
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turbulent  shear  stress  equation  reduces  to  the  form,  again  in  the  nomenclature  of  the 
present  analysis, 


r  turb 


(29) 


which  is  identical  with  the  results  of  Eqs.  (25)  and  (26)  derived  above  for  the  turbulent 
shear  contributions.  Arguments  for  the  preference  of  either  variance  or  invariance  of  the 
turbulent  shear  stress  vector  relative  to  the  mean  rate  of  strain  direction  are  currently 
based  on  slender  evidence.  The  analytical  studies  by  Mellor  (Ref.  27)  and  Cooper  (Ref. 
34)  concerning  three-dimensional  incompressible  turbulent  boundary  layers  indicate  that 
the  above-defined  invariant  model  of  turbulence  results  in  very  acceptable  calculated 
boundary-layer  parameters  relative  to  experiment.  Much  more  work,  especially  carefully 
controlled  experiments  involving  three-dimensional  compressible  turbulent  boundary  layers, 
remains  to  be  done  before  the  question  of  variance  or  invariance  can  be  completely  resolved. 
Until  such  time,  the  above  invariant  turbulence  model  appears  plausible  for  studying 
three-dimensional  turbulent  boundary-layer  flow. 

2.4  MIXING  LENGTH  MODEL 

The  turbulent  shear  stress  in  a  three-dimensional  turbulent  boundary  layer  as  governed 
by  Eqs.  (25)  and  (26)  is  treated  herein  by  the  use  of  a  two-layer  inner-outer  model  using 
Prandtl's  mixing  length  hypothesis  and  a  modification  of  van  Driest's  analysis  for  the  near 
wall  region.  This  results  in  a  continuous  distribution  of  the  shear  stress  from  the  laminar 
value  at  the  wall,  through  the  fully  turbulent  region,  reaching  zero  at  the  outer  edge 
of  the  boundary  layer.  The  energy  transport  in  a  turbulent  boundary  layer  is  treated  in 
this  work  through  the  incorporation  of  the  eddy  conductivity,  k,  into  a  turbulent  Prandtl 
number,  Prt.  Exactly  the  same  empirical  functions  are  used  as  in  the  two-dimensional 
flow  reported  by  Adams  (Ref.  11).  This  is  justified  by  noting  that  the  scalar  properties 
of  a  turbulence  field  are  unlikely  to  be  affected  by  moderate  three-dimensionality  because 
turbulence  is  inherently  three-dimensional  in  nature  for  even  so-called  two-dimensional 
flows. 

After  Escudier  (Ref.  35),  Patankar  and  Spalding  (Ref.  3)  recommend  the  following 
variation  of  the  mixing  length,  8*,  across  the  turbulent  two-dimensional  boundary  layer 
which  is  adopted  for  the  three-dimensional  case  per  the  above  discussion: 

=  k*y,  forO  <  y  =  Ayg/k* 

=  Ayg,  f°r  Ayg/k,  <  y 

where  the  values  for  the  various  numerical  constants  are  taken  to  be  k*  =  0.435  and 
X  =  0.09.  The  value  of  y  at  the  point  where  the  velocity  in  the  boundary  layer  is  equal 
to  0.99  of  the  velocity  at  the  boundary  layer  outer  edge  is  used  to  define  the  distance 
y£.  The  above  choices  follow  Patankar  and  Spalding  and  result  in  good  agreement  with 
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the  mixing  length  model  proposed  by  Maise  and  McDonald  (Ref.  36)  for  compressible 
two-dimensional  turbulent  boundary  layers. 

By  analogy  with  Stokes'  solution  for  an  infinite  flat  plate  undergoing  simple  harmonic 
motion  parallel  to  itself  in  an  infinite  fluid,  van  Driest  (Ref.  37)  concluded  that  in  the 
vicinity  of  a  wall  the  total  shear  stress  in  a  turbulent  two-dimensional  fluid  should  be 
of  the  form 


T 


+  P 


k*V 


-  exp 


(31) 


which  results  in  an  exponential  damping  of  the  turbulent  part  of  the  shear  stress  as  the 
wall  is  approached  and  yields  exactly  the  laminar  shear  stress  form,  r  =  pt(3u/dy),  at  the 
wall.  Although  Eq.  (31)  was  originally  developed  for  incompressible  flow,  it  can  be  applied 
to  compressible  flow  by  application  of  the  suggestion  by  Patankar  and  Spalding  (Ref. 
3)  that  the  local  value  of  shear  stress  be  used  instead  of  the  wall  value  as  originally 
recommended  by  van  Driest  (Ref.  37).  Hence,  by  analogy  of  Eq.  (31)  with  Eqs.  (25) 
and  (26),  the  relationship  for  the  three-dimensional  near-wall  shear  stress  as  used  in  the 
present  analysis  is 


T 


(32) 


where  the  constant  A,  is  taken  to  be  26.0  following  the  original  van  Driest  proposal 
(Ref.  37).  Note  that  the  damping  term  in  Eq.  (32)  reflects  the  application  of  the  local 
total  shear  stress  as  opposed  to  the  wall  shear  stress  of  Eq.  (31)  as  discussed  previously. 

Based  on  Eqs.  (26),  (27),  (30),  and  (32),  the  eddy  viscosity  expression  for  the  inner 
region  is 


while  for  the  outer  region 


_  .  2  2  <5G 

=  P  A2y£z  ^  (34) 

while  the  constants  k* ,  A* ,  A,  and  yfi  as  defined  previously.  The  constraint  used  to  define 
the  end  of  the  inner  region  and  the  beginning  of  the  outer  region  is  the  continuity  of 
the  eddy  viscosity.  From  the  wall  outward,  the  expression  for  the  inner  eddy  viscosity 
applies  until  ej  =  e0  from  which  point  the  outer  eddy  viscosity  is  used.  A  schematic 
of  this  variation  in  terms  of  the  mixing  lengths  is  shown  below. 
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Constant 


The  turbulent  Prandtl  number  defined  as 


(35) 


is  physically  a  measure  of  the  ratio  of  eddy  viscosity  to  eddy  conductivity,  that  is,  the 
ratio  of  the  turbulent  transport  of  momentum  to  the  turbulent  transport  of  heat.  Since 
the  flow  in  the  outer  region  of  a  turbulent  boundary  layer  shows  some  similarity  to  a 
turbulent  wake  flow,  one  may  advance  arguments  that  a  realistic  formulation  of  turbulent 
Prandtl  number  requires  a  separate  expression  for  the  inner  and  outer  region  just  as  in 
the  eddy  viscosity  formalism.  Recent  reviews  by  Rotta  (Refs.  38  and  39)  and  Meier  and 
Rotta  (Ref.  40)  indicate  that  the  turbulent  Prandtl  number  varies  in  magnitude  from 
approximately  0.8  near  the  boundary-layer  outer  edge  to  approximately  1.5  near  the  wall. 
Much  experimental  work  remains  to  be  done  in  defining  the  turbulent  Prandtl  number 
distribution  as  the  reviews  by  Rotta  point  out.  In  the  present  work  the  turbulent  Prandtl 
number  defined  by  Eq.  (35)  is  taken  to  remain  constant  at  the  value  0.90  across  the 
entire  boundary  layer  as  recommended  by  Patankar  and  Spalding  (Ref.  3)  for 
two-dimensional  turbulent  boundary  layers.  This  choice  of  turbulent  Prandtl  number  is 
consistent  with  the  analyses  of  Patankar  and  Spalding  (Ref.  3),  Martellucci,  Rie,  and 
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Sontowski  (Ref.  6),  Mayne  and  Dyer  (Ref.  7),  and  Adams  (Ref.  11).  The  early  analysis 
by  Smith  and  Cebeci  (Ref.  2)  assumed  a  value  of  unity,  whereas  the  recent  work  by 
Cebeci  (Ref.  8)  employs  a  variable  turbulent  Prandtl  number  model.  Bushnell  and  Beckwith 
(Ref.  4),  as  well  as  Hunt,  Bushnell,  and  Beckwith  (Ref.  29),  also  use  a  variable  turbulent 
Prandtl  number  model. 

2.5  GOVERNING  BOUNDARY-LAYER  EQUATIONS  IN  PLANE  OF  SYMMETRY 


In  order  to  initiate  numerical  solution  of  the  complete  three-dimensional 
boundary-layer  equations  using  a  marching,  implicit  finite-difference  procedure  to  integrate 
the  parabolic  partial  differential  equations,  initial  boundary-layer  profiles  are  required  along 
an  initial  or  starting  boundary.  For  the  case  of  a  body  of  revolution  at  angle  of  attack, 
the  initial  or  starting  boundary  is  the  windward  plane  of  symmetry.  For  the  sharp  cone 
geometry  of  interest  in  the  present  work,  the  windward  plane  of  symmetry  corresponds 
to  the  most  windward  ray  of  the  sharp  cone;  i.e.,  <j>  =  0  in  Fig.  1.  At  this  symmetry 
plane,  the  crossflow  velocity  component,  w,  vanishes  due  to  symmetry.  In  addition,  the 
crossflow  pressure  gradient,  9p/30,  also  vanishes  due  to  symmetry.  Under  these  restrictions, 
and  under  the  eddy  viscosity  model  adopted  in  Section  2.3,  every  term  in  the 
circumferential  momentum  equation  (3)  vanishes  identically.  Noting  that  under  these 
restrictions  the  continuity  equation  (1)  still  contains  the  crossflow  velocity  gradient,  9w/90, 
which  does  not  vanish  in  the  plane  of  symmetry,  Moore  (Ref.  16)  proposed  to  differentiate 
the  crossflow  equation  (3)  with  respect  to  <j>  first,  and  only  then  to  neglect  the  terms 
which  vanish  in  the  plane  of  symmetry.  This  procedure  applied  to  Eq.  (3)  yields  the 
following  results,  evaluated  in  the  plane  of  symmetry: 

CIRCUMFERENTIAL  (<f>)  MOMENTUM 


where 


+  p  V 


<5*0 

w 


+ 


Pu 


(36) 


_  dv/ 


(37) 


In  the  plane  of  symmetry  where  w  =  0,  the  continuity,  x-momentum,  and  energy  equations 
(1),  (2),  and  (9),  respectively,  reduce  to  the  following  form: 

CONTINUITY 


■k(pur)  +  £(pVr) +  pw0  = 0 

STREAMWISE  (x)  MOMENTUM 


_ dn 

pndZ 


+ 


+ 


(38) 


(39) 
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ENERGY 


put,  + 


,<5h 

"Vy- 


- 

u  si 


dp  /  v/du  V  d  7  /i  (  \  dhl 

d"x  +  (/1  +  f\dV  J  +  d^  [l,P7  +  PrJ  dyj 


(40) 


The  eddy  viscosity  in  the  above  equations  is  evaluated  from  Eq.  (27)  as 


— p2  I  tfu  I 

'  "  pl-  !3Fl 


(41) 


in  the  plane  of  symmetry.  The  absolute  magnitude  bars  on  du/dy  in  Eq.  (41)  take  into 
account  that  the  sign  on  rturb  must  change  with  du/dy  as  in  Eq.  (21)  presented  previously. 

At  this  point  it  is  well  to  note  that  the  plane  of  symmetry  restriction  has  reduced 
the  x-momentum  and  energy  equations  to  a  quasi-two-dimensional  form  in  that  they  do 
not  explicitly  depend  on  the  crossflow  velocity  gradient,  W0.  The  continuity  equation 
(38)  serves  to  couple  implicitly  the  streamwise  momentum  and  energy  equations  with 
the  crossflow  momentum  equation.  The  crossflow  velocity  gradient  term,  pw^,  in  the 
continuity  equation  (38)  can  be  physically  considered  as  a  local  source  (outflow)  or  local 
sink  (inflow)  term  with  the  strength  and  development  history  of  the  source  or  sink 
controlled  by  the  crossflow  momentum  equation  (36).  This  concept  is  especially  helpful 
in  understanding  how  a  quasi-two-dimensional  flow  reflects  three-dimensional  influence. 

Equation  (41)  shows  that  the  plane  of  symmetry  restriction  has  resulted  in  the  eddy 
viscosity  relationship  becoming  quasi-two-dimensional  in  form.  The  crossflow  influence  is 
reflected  implicitly  through  the  mixing  length,  £*,  which  is,  in  turn,  solely  dependent 
on  local  boundary-layer  conditions.  In  this  respect  the  application  of  the  local  total  shear 
stress  in  Eq.  (33)  is  especially  important. 

2.6  COORDINATE  TRANSFORMATION 

In  order  to  facilitate  numerical  integration  of  the  governing  quasi-two-dimensional 
plane  of  symmetry  boundary-layer  equations  (36),  (38),  (39),  and  (40),  it  is  convenient 
to  transform  them  to  a  coordinate  system  that  removes  the  singularity  at  x  =  0  and  stretches 
the  coordinate  normal  to  the  surface,  as  is  usually  done  in  two-dimensional  laminar  flow. 
The  coordinate  transformation  used  in  the  present  work  is  the  well-known 
Lees-Dorodnitsyn  (Ref.  41,  p.  31)  transformation.  The  new  independent  variables 
introduced  are 


£(*)  =  /  x  pefJ.  eUer2dx 
0 


(42) 


i/x.)) 


^y^dv 

\/2£  Jo  Pe 


(43) 
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evaluated  in  the  plane  of  symmetry.  The  transformed  streamwise  and  normal  derivatives 
hence  become 


d  TI  2  <5  drj  <9 

57  =  pe^euer  5?  +  57  5^ 

(44) 

(9  P U  er  <9 

dY  ~  s/2£  dr> 

(45) 

Introducing  the  two  scalar  stream  functions,  ^(x,y)  and  i^(x,y),  in  such  a  manner  as  to 
identically  satisfy  the  continuity  equation  (38)  following  the  approach  by  Moore  (Ref. 
13), 

dV 

pur  =  -t— 

ey  (46) 

difi 

p**1  =  57 

(47) 

(48) 

and  defining  transformed  stream  functions  f({,tj)  and  g(£,i?)  such  that 

Y(x,y)  =  V5?  (49) 

^(x,y)  =  yfe  tf*)  (50) 

the  governing  plane  of  symmetry  boundary-layer  equations  (36),  (39),  and  (40)  become, 
in  the  transformed  (£,r?)  plane, 


STREAMWISE  (x)  MOMENTUM 

«*f")'  +  [f  +  Gg]f"  +  /3p-  -  (f')2]  =  f 

CIRCUMFERENTIAL  (0)  MOMENTUM 


(£  g") '  +  [f  *  Gg]  g  "  -  f'g'  £  +  G  -  G[(g')2  +  x$] 


(52) 


ENERGY 
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ft?  dTe1 

Te  "^J  (53) 
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with  the  new  dependent  variables 


ne.r,)  =  f- 


(54) 


=  -jr 

0  =  -rf- 

1  e 


(55) 

(56) 


The  following  definitions  apply  to  the  above  equations 


£* 

PePe'  P> 

(57) 

r 

-  p') 

PePe\  P  PV 

(58) 

ft 

dUe 

= 

(59) 

X 

1  <?2P 

ypMg  d<f>2 

(60) 

e 

dU 

-PU  — e 

Pe  e  dx 
dp 

(61) 

dx 

G 

(62) 

dx 

ue 

v/yRTe 

(63) 

where  the  inviscid  streamwise  plane  of  symmetry  momentum  equation 


dp 

dx 


peue 


(64) 


has  been  used  to  relate  the  inviscid  velocity,  Ue,  to  the  imposed  inviscid  static  pressure. 
For  the  classical  boundary-layer  analysis  where  the  inviscid  flow  body  streamline  results 
are  used  as  the  conditions  at  the  outer  edge  of  the  boundary  layer,  the  parameter  e  = 
1.  For  the  more  general  case  where  swallowing  of  the  inviscid  flow  into  the  boundary 
layer  is  taken  into  account,  as  must  be  done  for  blunt-nosed  bodies  in  hypersonic  flow, 
the  value  of  e  must  be  determined  according  to  Eq.  (61).  This  approach  is  necessitated 
because  Eq.  (64)  is  only  valid  along  an  inviscid  streamline.  The  use  of  the  e  parameter 
in  the  above  manner  follows  the  approach  by  Blottner  (Ref.  42). 
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The  parameter  x  defined  by  Eq.  (60)  above  may  be  physically  interpreted  as  a  driving 
potential  or  forcing  function  for  the  plane  of  symmetry  crossflow  momentum  equation 
(52).  It  is  controlled  solely  by  the  inviscid  flow  and  must  be  input  to  the  analysis  along 
with  the  other  outer  edge  conditions  for  the  boundary  layer.  Furthermore,  since  x  is 
directly  dependent  on  the  second  derivative  of  the  pressure  distribution  in  the  ^-direction, 
information  concerning  the  allowable  boundary-layer  development  proceeding  away  from 
the  plane  of  symmetry  is  communicated  through  this  term.  In  short,  x  is  the  key  and 
controlling  parameter  in  determining  the  crossflow  influence  through  the  crossflow  velocity 
gradient  w<£  in  a  plane  of  symmetry  analysis. 


The  physical  boundary  conditions  given  by  Eqs.  (10)  and  (11)  become,  in  terms 
of  the  transformed  variables  in  the  plane  of  symmetry, 


MOMENTUM 


f(£if  =  0) 
f'(^  =  0) 
lim 

Tf->oo 

lim 


g  (£»)7  =  0)  =  0 
=  0)  =  0 


g '(£*?) 


=  l 


II 


=  g. 


(65) 


ENERGY 


0(£,7]  =  0) 


lim 

Tf-Mxj 


(66) 


The  outer  edge  boundary  condition,  g'e ,  for  the  crossflow  momentum  equation  can 
be  derived  as  follows.  From  the  asymptotic  character  of  the  circumferential  momentum 
equation  (52)  as  the  outer  edge  is  approached  from  below 

lim  g  "<£>?)  =  o 

7}-*oo 

lim  ,,,*  >  \  « 

g  (£»?)  =  o 

Tf-teo 


and 


lim 

7]-too 


=  g'e 
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so  that,  with  the  use  of  the  boundary  conditions  (65)  and  (66),  Eq.  (52)  reduces  to 
the  algebraic  form 


<<o2  +  |?+GSe'.  * G*  ■ 0  <«7> 

Solving  this  quadratic  equation  for  g'e  yields 

4 + gtx)  *  4+g£)!-  4°2x  <«*> 

«e  =  - 2G - 

The  present  study  is  concerned  solely  with  analysis  of  the  windward  plane  of  symmetry 
along  a  sharp  cone  at  angle  of  attack.  For  this  flow  condition  the  appropriate  choice 
of  sign  in  Eq.  (68)  is  positive  (+)  which  can  be  easily  determined  by  noting  that  g'e 
must  be  zero  for  the  sharp  cone  at  zero  angle  of  attack  where  x  is  also  zero  due  to 
the  uniform  circumferential  pressure  distribution.  Along  the  windward  plane  of  symmetry 
for  the  sharp  cone  at  angle  of  attack,  x  is  always  a  negative  quantity  so  that  g'c  is  always 
positive  denoting  a  source  (outflow)  influence  on  the  governing  equations. 

2.7  TRANSITION  ZONE  DESCRIPTION 

The  distance  between  transition  onset  at  x  =  Xt  and  the  beginning  of  fully  turbulent 
flow  further  downstream  at  x  =  Xt  is  denoted  herein  as  the  transition  zone.  In  the  present 
work,  which  is  restricted  solely  to  analysis  of  sharp  cone  windward  streamline  flows,  no 
attempt  is  made  to  predict  or  correlate  a  priori  the  location  of  transition  onset  at  Xt. 
Instead,  the  experimentally  indicated  beginning  of  the  transition  region  is  used  for  the 
onset  point  Xt  based  on  the  particular  sharp  cone  data  under  examination. 

Repeatability  of  gross  characteristics  of  the  transition  region  such  as  a  nondimerisional 
length  has  been  observed  in  experimental  investigations.  Potter  and  Whitfield  (Ref.  43) 
show  that  the  transition  zone  local  Reynolds  number,  Rec  Ax  (where  Ax  =  Xt  -  Xt), 
has  a  dependence  on  both  transition  onset  local  Reynolds  number,  ReC|xt,  and  the  local 
Mach  number,  Me.  In  the  present  work,  the  extent  of  the  transition  zone  will  be  represented 
by  a  transition  zone  local  Reynolds  number  ratio,  ReCjxx/Reeixt>  which  is  related  to 
Ree,Ax  by 


Re 


e,  Xr 


Re 


=  1 


5,  A: 


e,  X 


Re 


e,  X. 


(69) 


where  Reeixx  is  the  local  Reynolds  number  at  the  inception  of  fully  turbulent  flow. 
The  onset  of  transition  in  the  present  work  is  determined  from  examination  of  experimental 
heat-transfer  data.  The  end  of  transition  is  taken  herein  to  coincide  with  the  location 
of  experimentally  determined  maximum  heating.  A  schematic  of  this  nomenclature  is  shown 
on  the  following  page. 
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Heat 

Transfer 


Surface  Distance  along  Body 


Masaki  and  Yakura  (Ref.  12)  present  a  comprehensive  survey  of  unclassified 
experimental  data  showing  transition  on  flat  plates  and  cones  over  a  range  of  local  Mach 
numbers  from  2  to  15.  Based  on  this  survey  (which  is  given  in  Fig.  2  of  their  paper) 
they  suggest  that  a  value  for  Ree,xT/Ree,x t  of  approximately  two  appears  reasonable 
at  high  supersonic  speeds.  The  value  of  two  is  not  a  universal  number,  and  a  great  deal 
of  work  remains  to  be  done  before  this  boundary-layer  characteristic  is  fully  understood. 
The  effects  of  wall  cooling,  pressure  gradient,  and  three-dimensional  crossflow  are  moot 
points  at  this  time.  It  should  be  noted,  however,  that  Murphy  and  Rubesin  (Ref.  44) 
found  that  Re^xj/Ree.Xt  *  2  based  on  flight  test  results  from  the  Mark-2  blunt-nosed 
ballistic  reentry  vehicle.  In  addition,  the  recent  experimental  data  reported  by  Young, 
Reda,  and  Roberge  (Ref.  45)  taken  in  a  hypersonic  wind  tunnel  (Mach  10)  on  a  lifting 
entry  configuration  indicate  that  the  factor  of  two  ratio  of  transition  end-to-onset 
boundary-layer-edge  Reynolds  number  provides  a  good  description  of  transition  zone  extent 
for  even  complex  geometries  at  angle  of  attack.  As  the  transition  onset  Reynolds  number 
is  increased,  the  ratio  of  transition  end-to-onset  boundary-layer-edge  Reynolds  number 
appears  to  slightly  decrease  below  the  value  of  two  (see  Refs.  12,  45,  and  46  for  further 
discussion  of  this  point). 
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Based  on  the  above  discussion,  the  present  investigation  utilizes  a  transition  zone 
description  following  the  criterion 


Re 


e,  X. 


Re 


=  2 


e,  X. 


(70) 


which  is  identical  to  that  used  previously  by  Adams  (Ref.  1 1 )  with  success. 


2.8  INTERMITTENCY  FACTOR  TREATMENT  OF  TRANSITION 


In  order  to  provide  a  suitable  progressive  transition  model  as  opposed  to  an 
instantaneous  transition  model,  the  present  analysis  utilizes  the  phenomenological  theory 
of  the  transition  zone  developed  by  Emmons  (Refs.  47  and  48)  which  is  based  upon 
the  experimental  observation  that  transition  in  a  boundary  layer  is  characterized  by  the 
intermittent  appearance  of  turbulent  spots  which  grow  and  move  downstream  with  the 
fluid.  In  the  transition  zone,  the  intermittent  appearance  of  turbulent  spots  and  their 
growth  as  they  move  downstream  produces  alterations  of  laminar  and  turbulent  flow  with 
a  gradual  increase  in  the  mean  turbulence  level.  At  any  point  in  the  transition  zone,  the 
probability  of  turbulent  flow  is  prescribed  by  a  model  based  on  the  residence  times  of 
intermittent  turbulent  spots.  Dhawan  and  Narasimha  (Ref.  49)  modified  Emmons'  original 
formulation  of  the  probability  of  turbulent  flow  to  account  for  a  localized  region  of 
breakdown  about  x  =  Xt  in  conjunction  with  a  single  universal  intermittency  distribution. 
The  Dhawan  and  Narasimha  model  simply  expresses  the  intermittency  factor.  If,  which 
indicates  the  probability  of  the  flow  being  turbulent  at  some  point  as 


where  $  is  a  constant  equal  to  0.412  and  x  is  a  measure  of  the  extent  of  the  transition 
region  characterized  by 


x  =  =  0.75”  ^lf  =  0.25  (72) 

The  beginning  of  the  transition  zone,  where  If (0)  =  0,  was  established  by  Dhawan 
and  Narasimha  (Ref.  49)  as  being  located  at  the  transition  onset  distance  x  =  Xt.  In 
order  to  describe  the  extent  of  the  transition  zone,  Masaki  and  Yakura  (Ref.  1 2)  assume 
that  fully  turbulent  flow  is  achieved  at  x  =  Xj  where 

If(XT)  £  0.97  (73) 

Substituting  this  assumption  into  Eq.  (71)  gives  the  length  of  the  transition  zone  as 

Ax  =  XT  -  Xt  =  2.96 x  (74) 

In  order  to  evaluate  x>  recall  from  Eq.  (70)  that 

Ree,  XT  =  2Ree,  xt  (75) 


19 


AEDC-TR-71-235 


which  implies  physically  that  the  transition  zone  length  Ax  is  equal  to  the  transition  onset 
distance  Xt.  Hence  this  reveals  that 

Ax  =  xt  (76) 


so  that  x  becomes,  from  Eq.  (74), 


X. 

*  =  2M 


(77) 


Under  these  restrictions  the  intermittency  factor  given  by  Eq.  (71)  takes  the  final  form 


I,(x)  =  1 


exp 


-0.412  (2.96)2(j-  -  l) 


where  x  =  Xt 


(78) 


Although  the  original  derivation  of  the  intermittency  factor  was  based  on  experiments 
performed  in  low  speed  incompressible  flow,  Dhawan  and  Narasimha  (Ref.  49)  established 
its  validity  for  use  in  prediction  of  supersonic  skin  friction  on  a  flat  plate.  To  the  author's 
knowledge,  no  work  has  been  done  on  the  effects  of  pressure  gradients  or  crossflow  on 
the  If(x)  distributions  in  either  incompressible  or  compressible  flows,  although  Dhawan 
and  Narasimha  speculate  that  pressure  gradients  would  chiefly  influence  the  manner  of 
spot  growth  in  the  initial  period  near  If  =  0.  Another  point  of  interest  is  the  distribution 
of  If  across  the  boundary  layer,  viz.,  normal  to  the  surface.  Dhawan  and  Narasimha  state 
that,  although  this  If(y)  variation  is  probably  of  importance  to  the  detailed  structure  of 
the  turbulent  motion,  the  value  of  If  near  the  wall  is  the  characteristic  property  of 
importance  for  the  transition  region  and  the  If(y)  variation  has  only  a  secondary  influence 
in  determining  the  mean  velocity  profiles  in  transition.  Hence,  for  the  supersonic  and 
hypersonic  sharp  cone  windward  streamline  flows  of  interest  in  the  present  investigation, 
Eq.  (78)  is  used  for  the  intermittency  factor  distribution  in  order  to  assess  its  validity 
for  use  in  such  regimes. 


In  the  present  work  the  intermittency  factor  multiplies  the  eddy  viscosity  in  the 
modified  density-viscosity  product  ratios  £*  and  £**  defined  by  Eqs.  (57)  and  (58)  in 
the  following  manner: 


£* 


(79) 


r  _  _ej£-(i  +  it'  £:) 


(80) 


The  physical  justification  for  using  the  intermittency  factor  in  this  manner  lies  in  the 
basic  physical  definition  of  the  intermittency  factor  as  representing  the  fraction  of  time 
any  point  spends  in  turbulent  flow,  i.e.,  the  probability  of  the  flow  being  turbulent  at 
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a  given  point.  If  If  =  0,  the  flow  is  completely  laminar,  whereas  if  If  =  1 ,  the  flow  is 
fully  turbulent.  The  present  treatment  of  the  intermittency  factor-eddy  viscosity  product 
in  the  governing  equations  simply  means  that  the  (fully  turbulent)  eddy  viscosity  is 
multiplied  by  a  damping  factor  (the  intermittency  factor)  in  order  to  achieve  a  smooth 
and  continuous  variation  from  laminar  to  turbulent  flow  through  the  transition  region. 
Similar  approaches  have  been  reported  by  Adams  (Ref.  11),  Harris  (Ref.  10),  and 
Martellucci,  Rie,  and  Sontowski  (Ref.  6).  These  references  indicate  that  use  of  the 
intermittency  factor  in  the  above-defined  manner  results  in  very  acceptable  calculation 
of  transitional  heating  rates  when  applied  in  conjunction  with  the  factor  of  two  transition 
zone  extent  criterion  discussed  previously.  Certainly  more  work  remains  to  be  done  relative 
to  the  definition  and  application  of  the  intermittency  factor  concept  in  three-dimensional 
flows. 

2.9  BOUNDARY-LAYER  PARAMETERS 

Given  the  numerical  solution  to  the  governing  equations  of  motion  (Eqs.  (51),  (52), 
and  (53))  following  the  integration  procedure  of  Appendix  III,  the  associated  local 
boundary-layer  parameters  at  a  given  body  station  may  be  determined  as  follows.  The 
local  convective  heat  flux  at  the  body  surface  (y  =  0)  is  given  by  the  well-known  Fourier 
law, 

^  =  _k"(|y')V¥  (81) 


which  becomes,  in  the  transformed  (£,i?)  coordinates, 


-We^pT, 

PrwV^ 


0'(£> 7  =  0) 


(82) 


where 


n  _  W  (83) 

W  Pe^e 

In  a  similar  manner,  the  local  wall  shearing  stress  along  the  windward  streamline  may 
be  written  in  terms  of  the  transformed  coordinates  as 


r 


w 


f”(£»?  =  0) 


(84) 


The  actual  height  of  the  boundary  layer  along  the  windward  streamline  is  determined 
from  the  y-transformation  relationship,  Eq.  (43),  as 


Y 


PeUe  r 


(85) 
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The  classical  definition  of  the  two-dimensional  boundary-layer  displacement  thickness  is 
given  by 


(86) 


which  becomes,  in  transformed  coordinates, 

s*  =  f-ie-nir,  (87) 

PeUer 

For  the  case  of  boundary-layer  flow  along  the  windward  streamline  of  a  sharp  cone 
at  incidence,  two  lengths  characterizing  mass-flow  defect  may  be  defined  in  terms  of  the 
profiles  of  the  two  velocity  components  tangential  to  the  surface  (see  Fig.  1  for  clarification 
of  nomenclature): 


and 


(88) 


(89) 


and  is  is  not  clear  which,  if  either,  defines  a  displacement  thickness  properly  describing 
the  extent  to  which  the  nonviscous  flow  is  deflected  by  the  boundary  layer.  Moore  (Ref. 
50)  has  considered  this  problem  in  great  detail  relative  to  the  proper  definition  of  a 
three-dimensional  displacement  surface.  For  the  special  case  of  windward  streamline 
boundary-layer  flow  on  a  sharp  cone  at  incidence,  Moore's  Eq.  (14)  on  page  7  of  Ref. 
50  shows  that 


A 


* 


(90) 


where  A*  is  the  three-dimensional  plane  of  symmetry  boundary-layer  displacement 
thickness.  Note  that  A*  differs  from  5X*  in  the  plane  of  symmetry  due  to  the 
circumferential  mass-flow  defect.  In  terms  of  the  transformed  (|,7j)  coordinates,  the  two 
mass-flow  defect  lengths  may  be  written  as 


S*  =  'fit  (0-  i')drj 

X  PetjerJo 


(91) 


and 


valid  along  the  windward  streamline. 


(92) 
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Various  nondimensional  quantities  of  interest  are  defined  as  follows: 


LOCAL  REYNOLDS  NUMBER  BASED  ON  FREE-STREAM  CONDITIONS 


Re. 


p  U  X 
r  oo  sc 


OO,  X 


(93) 


LOCAL  REYNOLDS  NUMBER  BASED  ON  INVISCID  EDGE  CONDITIONS 


Re 


e,  x 


(94) 


LOCAL  SKIN-FRICTION 
CONDITIONS 


COEFFICIENT  REFERENCED  TO  FREE-STREAM 


C,  = 


-p  U 2 

2  r  oo  oo 


(95) 


LOCAL  STANTON  NUMBER  BASED  ON  INVISCID  EDGE  CONDITIONS  AND 
ADIABATIC  WALL  TEMPERATURE 


St. 


PeUeCp(Taw-Tw) 


(96) 


LOCAL  STANTON  NUMBER  BASED  ON  FREE-STREAM  CONDITIONS  AND 
ADIABATIC  WALL  TEMPERATURE 


St 


-flw 

p  u  C  (T  TTl 
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(97) 


In  the  above 


T 

=  r{  +  7= -  (1  -  rf) 

1  o,  oo 


(98) 


with  rf  the  recovery  factor  which  must  be  defined  relative  to  the  flow  under  examination, 
i.e.,  for  laminar  boundary  layers 


rf  =  Vp7  =  v'frJI  =  0.843  (99) 

and  for  turbulent  boundary  layers 

rf  =  Vp7  =  3N/oTl  =  0.892  (100) 

where  Pr  =  0.71,  as  previously  specified  in  Section  2.1. 
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2.10  NUMERICAL  SOLUTION  OF  THE  GOVERNING  BOUNDARY-LAYER 
EQUATIONS 

In  the  present  work  (restricted  to  the  windward  plane  of  symmetry  of  a  sharp  cone 
at  incidence),  numerical  solution  of  the  simultaneous,  nonlinear,  parabolic,  partial 
differential  equations  (Eqs.  (51),  (52),  and  (53))  is  performed  by  obtaining  linear 
finite-difference  equivalents  of  the  equations  and  solving  these  using  an  iterative,  marching, 
implicit  finite-difference  technique  involving  inversion  of  tridiagonal  matrices.  Full  details 
of  this  numerical  approach  are  given  in  Appendix  III. 

The  digital  computer  code  is  written  in  FORTRAN  63  for  use  on  a  CDC  1604-B 
digital  computer.  Solution  time  including  printout  is  approximately  20  minutes  total  on 
the  1604-B  machine  for  a  windward  streamline  divided  into  100  stations.  No  numerical 
stability  problems  have  been  encountered  with  the  present  finite-difference  approach  due 
to  its  implicit  nature. 

2.11  INVISCID  FLOW  ABOUT  A  SHARP  CONE  AT  INCIDENCE 

A  recent  investigation  by  Jones  (Ref.  51)  has  resulted  in  an  accurate  and  efficient 
numerical  integration  procedure  for  solution  of  the  governing  partial  differential  equations 
describing  the  supersonic  or  hypersonic  inviscid  flow  field  around  a  sharp  cone  at  incidence. 
Basically,  Jones'  method  uses  the  condition  of  conicity  to  reduce  the  problem  to  a  set 
of  elliptic  nonlinear  partial  differential  equations  in  two  independent  variables.  A 
transformation  of  coordinates  is  used  to  fix  the  boundaries,  one  of  which  is  the  unknown 
shock  wave,  between  which  the  elliptic  equations  are  to  be  satisfied.  This  transformation 
also  has  the  effect  of  including  the  body  shape  in  the  coefficients  of  the  partial  differential 
equations  and  in  the  boundary  conditions,  so  that  the  same  method  can  be  used  for  general 
conical  body  shapes  simply  by  changing  a  few  program  statements  to  redefine  the  equation 
of  the  body.  In  fact  the  method  is,  in  many  cases,  only  limited  by  locally  supersonic 
crossflow  conditions,  by  the  entropy  singularity  moving  too  far  away  from  the  surface, 
or  by  the  shock  approaching  very  close  to  the  Mach  wave.  In  practice,  these  restrictions 
limit  the  allowable  angle-of-attack  range  to  a/Sv  ^  1.2  (see  Fig.  1  for  clarification  of 
nomenclature). 

At  the  present  time  the  method  has  been  used  successfully  for  circular  cones  and 
for  bodies  that  can  be  obtained  by  successive  perturbations  of  a  circular  cone  and  that 
do  not  have  curvatures  that  are  too  large.  Jones  (Ref.  51)  has  reported  examples  for 
circular  cones  at  incidence,  elliptic  cones,  and  a  body  whose  cross-sectional  shape  is 
represented  by  a  fourth  order  even  cosine  Fourier  series. 

The  method  is  efficient  in  computer  time  compared  with  other  fully  numerical 
techniques,  and  one  solution  takes  from  about  one-half  minute  to  three  minutes  on  an 
IBM  360/50  computer  for  the  circular  cone  at  incidence  -  the  time  increasing  as  the 
incidence  increases.  This  is  to  be  compared  with  a  time  requirement  of  approximately 
one-half  hour  on  an  IBM  360/50  computer  for  the  technique  developed  by  Moretti  (Ref. 
52)  in  which  the  flow-field  solution  is  obtained  by  marching  step  by  step  downstream 
(approximately  400  downstream  steps  are  required)  until  a  conicity  condition  is  sufficiently 
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well  satisfied.  Comparison  of  results  between  the  Jones  and  Moretti  approaches  shows 
excellent  agreement,  with  the  Jones  digital  computer  code  being  a  factor  of  approximately 
ten  faster  than  the  Moretti  approach  in  solution  time.  An  analysis  very  similar  to  that 
of  Jones  has  recently  been  reported  by  South  and  Klunker  (Ref.  53),  and  Holt  and  Ndefo 
(Ref.  54)  have  developed  a  method  of  integral  relations  approach  to  the  problem.  The 
important  point  to  note  is  that  all  of  the  above-referenced  analyses  report  excellent 
agreement  with  experiment  for  sharp  circular  and  elliptic  cones  at  incidence  under 
supersonic  and  hypersonic  flow  conditions;  therefore,  the  choice  of  which  analysis  is  indeed 
the  best  remains  an  open  question.  The  present  author's  experience  with  use  of  the  Jones' 
digital  computer  code  (Ref.  55)  has  been  most  favorable  from  a  user’s  standpoint. 

For  the  present  investigation,  concerned  only  with  the  windward  plane  of  symmetry 
along  a  sharp  cone  at  incidence,  the  important  inviscid  parameters  required  for  input  to 
the  boundary-layer  analysis  as  outer-edge  conditions  are 

-  T 

— ,  ■=-  ,  M  ,  and  x 

p~  V  e’ 


where  recall  that 

_  _1 _  d2p 

ypMe  d4>2 
and 

g  e  2G 

from  Eqs.  (60)  and  (68).  respectively,  with  all  quantities  evaluated  on  the  windward 
streamline  where  0  =  0.  All  of  the  above-indicated  inviscid  parameters  may  be  determined 
from  the  Jones  digital  computer  code  with  numerical  differentiation  required  to  evaluate 
d-p/902  on  the  windward  streamline.  It  should  be  noted  that  the  so-called  "isentropic 
surface"  values  of  the  inviscid  parameters  (see  Ref.  51  for  clarification  of  this  terminology) 
have  been  used  in  the  present  investigation. 

It  should  be  pointed  out  in  conclusion  that  Jones  (Ref.  56)  has  recently  published 
a  very  complete  and  thorough  set  of  tables  for  inviscid  supersonic  and  hypersonic  flow 
about  circular  cones  at  incidence  in  a  perfect  gas,  y  =  1 .40,  stream.  These  tables  can 
be  used  to  provide  all  of  the  needed  inviscid  information  for  input  to  the  present 
boundary-layer  analysis. 


’SECTION  III 

RESULTS  AND  DISCUSSION 

As  pointed  out  by  Smith  (Ref.  57)  in  his  recent  review  of  boundary-layer  research 
over  the  past  decade,  the  ultimate  test  of  any  theoretical  analysis  of  the  turbulent  boundary 
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layer  is  the  end  accuracy  (and  not  the  intrinsic  logic)  since  empiricism  is  involved.  Hence, 
with  this  in  mind,  comparison  of  experimental  data  with  the  present  laminar,  transitional, 
and  turbulent  windward  plane  of  symmetry  boundary-layer  analysis  is  mandatory.  Since 
the  current  effort  is  limited  solely  to  sharp  cone  geometry,  the  following  additional 
restrictions  are  imposed  at  this  point  in  order  to  clearly  define  the  experimental  flow 
field  for  comparison  with  the  analytical  model: 

1.  Angle  of  incidence  range  0  —  a  %  1.2  6V. 

2.  Uniform  free  stream;  i.e.,  no  source  flow  effects. 

3.  Constant  uniform  wall  temperature. 

4.  Windward  streamline  transition  zone  well  defined  by  heat-transfer  data 
where  applicable. 

Under  the  above  restrictions,  a  review  of  unclassified,  published  literature  on  sharp  cone 
at  incidence  experimental  investigations  revealed  the  seven  sources  (Refs.  58  through  64) 
listed  in  Table  1,  Appendix  II.  Most  published  investigations  do  not  meet  the  rather  stringent 
restrictions  imposed  above  and  hence  must  be  eliminated  on  this  basis.  For  example,  Stetson 
and  Rushton  (Ref.  65)  present  no  usable  heat-transfer  data  defining  the  transition  zone 
at  angle  of  attack.  By  only  using  experimental  data  from  Table  I,  meaningful  comparisons 
between  theory  and  experiment  can  be  made  with  some  confidence  since  the  assumptions 
under  which  the  analytical  model  is  formulated  match  the  restrictions  placed  on  acceptable 
experimental  data. 

Table  II  lists  the  windward  streamline  inviscid  outer-edge  conditions  for  input  to 
the  boundary-layer  analysis  corresponding  to  the  various  flows  in  Table  I  for  which 
experimental  data  are  available.  The  inviscid  conditions  are  determined  from  numerical 
solution  of  the  inviscid  conical  flow  equations  for  a  sharp  cone  at  incidence  using  the 
Jones'  digital  computer  code  (Refs.  51  and  55)  described  in  Section  2.11.  Consistent  with 
the  gas  model  adopted  in  Section  2.1,  a  perfect  gas  with  y  =  1.40  is  used  in  the  inviscid 
analysis. 

In  order  to  assess  the  basic  applicability  of  the  present  windward  plane  of  symmetry 
boundary-layer  analysis  for  laminar  applications,  Figs.  2,  3,  4,  and  5  present  windward 
ray  heat-transfer-rate  calculations  using  the  present  theory  relative  to  experimental  data. 
Complete  flow  conditions  may  be  found  in  Table  I  with  the  corresponding  inviscid 
outer-edge  conditions  given  in  Table  II.  Figures  2  and  3  show  good  agreement  between 
theory  and  experiment  for  the  heat-transfer  ratio  qw(a)/qw  (a  =  0)  at  a  common  body 
location.  Figures  4  and  5  present  calculated  adiabatic  wall  Stanton  numbers  relative  to 
experimental  measurements.  Again,  the  agreement  between  theory  and  experiment  is  good 
except  for  the  front  half  of  the  distribution  presented  in  Fig.  5.  Included  in  Fig.  5  is 
a  zero-angle-of-attack  distribution  illustrating  the  results  from  an  intermittency  factor 
treatment  of  transition  using  the  present  theory  relative  to  experimental  data.  The  favorable 
heat-transfer  results  presented  above  indicate  that  indeed  the  present  theory  is  applicable 
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and  accurate  under  laminar  boundary-layer  conditions  along  the  windward  streamline  of 
a  sharp  cone.  The  case  of  transitional  and  turbulent  boundary  layers  along  the  windward 
plane  of  symmetry  will  now  be  considered. 

Of  the  data  in  Table  I,  the  investigation  by  McCauley,  Saydah,  and  Bueche  (Ref. 
62)  is  one  of  the  best  documented  with  respect  to  test  conditions  and  analysis  of  data. 
The  test  was  conducted  in  the  Hypersonic  Wind  Tunnel  (C)  of  the  von  KArmAn  Gas 
Dynamics  Facility,  Arnold  Engineering  Development  Center,  at  a  nominal  Mach  number 
of  10  under  cold  wall  conditions.  Complete  test  conditions  are  given  in  Table  I.  Figure 
6  shows  a  comparison  of  results  from  the  present  theory  relative  to  experimental 
heat-transfer  measurements  along  the  windward  streamline  under  natural  transition 
conditions.  In  general,  the  agreement  between  experimental  data  and  calculated  results 
from  the  present  theory  is  excellent  for  both  zero  and  nonzero  angles  of  incidence.  The 
transitional  heating  rates  are  well  defined  by  the  present  theory  which,  it  should  be  recalled, 
requires  that  the  location  of  onset  to  transition  be  prescribed  as  indicated  in  Fig.  6.  All 
three  sets  of  experimental  data  presented  in  Fig.  6  clearly  show  a  factor  of  about  two 
in  the  ratio  Xx/Xt  as  previously  discussed  in  Section  2.7  and  used  in  the  present  analytical 
model.  Note  from  Fig.  6  that  an  increase  in  angle  of  attack  caused  the  onset  of  transition 
location,  Xt,  to  move  slightly  rearward  along  the  windward  streamline  which  is  in  agreement 
with  the  findings  of  Stetson  and  Rush  ton  (Ref.  65).  Further  note  that  for  this  natural 
transition  condition  the  boundary-layer  flow  remains  transitional  over  the  entire  aft  half 
of  the  cone. 

The  primary  emphasis  in  the  experimental  investigation  of  Ref.  62  was  placed  on 
determining  the  effect  of  controlled  three-dimensional  roughness  (spherical  trips)  on 
hypersonic  laminar  boundary-layer  transition.  Presented  in  Fig.  7  is  a  comparison  of 
calculated  heat  transfer  relative  to  experimental  measurements  for  a  case  in  which  the 
laminar  boundary  layer  was  "tripped"  at  the  body  location  x/L  =  0.15.  For  this  condition 
the  calculations  assumed  an  instantaneous  transition  from  fully  laminar  to  fully  turbulent 
flow  at  the  above-specified  trip  location  on  the  body.  Also  shown  for  sake  of  comparison 
is  a  zero  angle-of-attack  calculation  using  the  intermittency  factor  treatment  of  transition 
with  onset  to  transition  taken  at  the  body  location  x/L  =  0.25.  The  important  point 
to  observe  from  Fig.  7  is  that  the  laminar  boundary  layer  is  not  effectively  "tripped" 
to  an  immediate  fully  turbulent  state  relative  to  the  numerical  calculations  for  the 
angle-of-attack  conditions.  Only  at  the  end  of  the  body  (x/L  =  1 .0)  has  the  experimental 
data  for  the  angle-of-attack  conditions  reached  the  level  predicted  by  the  present  theory. 
Note  especially  that  the  level  of  the  experimental  data  remains  essentially  constant  from 
the  trip  location  to  the  end  of  the  body  for  the  angle-of-attack  conditions.  For  zero  angle 
of  attack,  the  "tripped"  boundary-layer  response  is  sufficiently  rapid  to  produce  fully 
turbulent  flow  over  the  entire  aft  half  of  the  cone.  Compare  Figs.  6  and  7  for  clarification 
of  the  relative  behavior  between  "untripped"  and  "tripped"  boundary  layers. 

Turning  now  to  details  of  the  boundary-layer  structure  along  the  windward  plane 
of  symmetry,  Figs.  8  and  9  present  calculated  boundary-layer  profiles  (streamwise  velocity, 
crossflow  velocity,  static  temperature,  total  temperature,  Mach  number,  and  pitot  pressure) 
for  fully  laminar  (body  location  x/L  =  0.40)  and  fully  turbulent  (body  location  x/L  = 
1.0)  flow  conditions  corresponding  to  Fig.  6.  The  different  character  of  the  laminar  and 
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turbulent  profiles  is  apparent.  Note  the  overshoot  in  both  the  laminar  and  turbulent 
crossflow  velocity  profiles.  Direct  comparison  of  laminar  and  turbulent  profiles  at  a 
common  body  location  will  be  given  in  the  next  paragraph. 

In  order  to  present  the  differences  between  laminar  profiles  and  fully  turbulent  profiles 
at  a  given  common  body  location  along  the  windward  streamline,  Fig.  10  shows  calculated 
streamwise  velocity,  crossflow  velocity,  and  static-temperature  distributions  for  the  two 
types  of  flow  at  x/L  =  1.0.  Note  that  the  streamwise  velocity  profiles  for  turbulent  flow 
are  much  fuller  than  the  laminar  profiles,  as  expected.  Further  note  that  the  maximum 
overshoot  in  the  turbulent  crossflow  velocity  profile  is  less  than  the  corresponding 
angle-of-attack  laminar  value  which  means  physically  that  the  outflow  influence  along  the 
windward  streamline  of  a  laminar  boundary  layer  is  stronger  than  a  corresponding 
angle-of-attack  turbulent  boundary  layer.  The  implications  of  this  finding  as  it  affects 
surface  heat  transfer  will  be  discussed  in  the  next  paragraph.  With  respect  to  the 
static- temperature  distributions,  Fig.  10  shows  that  the  peak  static  temperature  in  the 
laminar  boundary  layer  occurs  much  further  from  the  wall  than  the  peak  turbulent  static 
temperature,  about  a  factor  of  three  in  distance  normal  to  the  surface.  This  factor  of 
three  in  location  of  peak  static  temperature  is  somewhat  reflected  in  the  thickness  of 
the  turbulent  boundary  layer  being  approximately  twice  the  thickness  of  the  laminar 
boundary  layer  as  can  be  seen  by  close  inspection  of  Fig.  10. 

A  concept  often  used  (see  Ref.  66,  for  example)  to  approximately  calculate  windward 
ray  heating  rates  under  angle-of-attack  conditions  is  the  so-called  "effective  cone"  approach 
in  which  a  zero  angle-of-attack  calculation  is  performed  on  an  "effective  cone"  which 
has  a  cone  half-angle  equal  to  the  physical  cone  half-angle  plus  the  physical  cone  angle 
of  attack.  Needless  to  say,  the  "effective  cone"  technique  does  not  properly  include  the 
effects  of  crossflow  (outflow)  as  it  affects  the  windward  ray  boundary-layer  structure. 
Presented  in  Fig.  1 1  arc  comparisons  of  "effective  cone"  calculations  relative  to  the  present 
three-dimensional  windward  ray  analysis  for  the  case  of  laminar,  transitional,  and  turbulent 
boundary^ layer  flow  at  various  angles  of  attack.  Note  that  the  flow  conditions  are  identical 
with  those  of  Fig.  6  discussed  previously.  Also,  onset  to  transition  locations  has  been 
matched  between  the  two  analyses  for  a  given  angle  of  attack.  As  can  be  seen  from  Fig. 
11,  both  the  heat  transfer  and  skin  friction  are  much  more  sensitive  to  three-dimensional 
crossflow  (outflow)  effects  for  a  laminar  boundary  layer  than  for  a  turbulent  boundary 
layer.  A  very  informative  summary  of  these  results  is  given  in  Fig.  12.  The  important 
point  here  is  that  "effective  cone"  calculations  appear  to  be  acceptable  for  the  present 
geometry  and  flow  conditions  provided  the  boundary  layer  is  fully  turbulent.  For  a  laminar 
boundary  layer,  "effective  cone"  calculations  result  in  a  20- to  30- percent  underprediction 
(based  on  the  x/L  =  0.40  laminar  results  presented  in  Fig.  12)  for  both  heat  transfer 
and  skin  friction  relative  to  the  present  three-dimensional  windward  streamline  analysis. 
These  results  confirm  the  statement  of  Vaglio-Laurin  (Ref.  24)  that  "due  to  the  larger 
shearing  stresses,  smaller  three-dimensional  effects  can  be  expected  for  turbulent  layers 
as  compared  with  laminar  layers  subject  to  the  same  boundary  conditions."  In  addition, 
the  recent  paper  by  Widhopf  (Ref.  66)  concludes  that  "equivalent  cone  techniques  do 
not  yield  uniformly  valid  approximations  but  may  provide  acceptable  results  in  some  cases." 
Indiscriminate  use  of  the  "effective  cone"  approach,  especially  for  laminar  boundary  layers, 
is  not  advised  by  the  present  author. 
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Details  of  the  boundary-layer  thickness  (based  on  the  normal  distance  from  the  surface 
where  f'  =  0.995)  and  the  three-dimensional  displacement  thickness  defined  by  Eq.  (90) 
are  presented  in  Fig.  1 3  for  the  same  flow  and  a ngle-of- attack  conditions  as  Fig.  6  previously 
discussed.  The  general  thinning  of  the  boundary  layer  as  the  angle  of  attack  increases 
is  apparent.  A  comparison  of  the  three-dimensional  displacement  thickness,  A*,  relative 
to  the  two-dimensional  streamwise  direction  value,  8X*,  defined  by  Eq.  (91)  is  given  in 
Fig.  14.  The  three-dimensional  displacement  thickness,  A*,  is  slightly  smaller  than  the 
two-dimensional  value  at  a  common  body  location.  Note  that  the  slopes  dA*/dx  and 
d8x*/dx  arc  very  similar  at  a  common  body  location  which  is  of  importance  in  computing 
the  displacement-induced  pressure  along  the  windward  streamline  using  an  inviscid 
tangent-cone  technique.  Further  note  that,  as  the  angle  of  attack  is  increased,  the  slope 
of  the  displacement  thickness  decreases  for  a  given  body  location.  This  means  physically 
that  displacement-induced  pressure  effects  will  become  of  less  and  less  significance  on 
the  windward  streamline  as  the  angle  of  attack  is  increased.  For  the  present  flow  conditions. 
Fig.  3  of  Ref.  62  shows  that  the  displacement-induced  pressure  is  negligible  even  at  zero 
angle  of  attack  and  hence  will  not  be  further  considered  in  the  present  work.  It  should 
be  pointed  out  in  passing,  however,  that  the  displacement-induced  pressure  increment  may 
attain  its  maximum  value  in  the  transition  region  as  can  be  seen  from  careful  examination 
of  the  displacement  thickness  slopes  obtained  from  Figs.  13  and  14. 

Turning  now  to  additional  comparisons  of  experimentally  obtained  transitional  and 
turbulent  boundary-layer  data.  Fig.  15  presents  a  comparison  of  a  heat-transfer  rate 
parameter  calculated  using  the  present  windward  plane  of  symmetry  theory  with  the  data 
reported  by  DiCristina  (Ref.  63)  for  Mach  10.2  hypersonic  flow.  This  test  was  conducted 
in  the  Hypersonic  Wind  Tunnel  (C)  of  the  von  K&rm&n  Gas  Dynamics  Facility,  Arnold 
Engineering  Development  Center;  test  conditions  are  given  in  Table  1  with  corresponding 
inviscid-edge  conditions  in  Table  II.  Agreement  of  results  from  the  present  theory  and 
experiment  is  excellent  with  the  transitional  zone  heating  well  defined  by  the  current 
intermittency  factor  analysis.  Note  the  change  in  recovery  factor  definition  between  laminar 
and  transitional  flow  which  is  consistent  with  the  manner  in  which  DiCristina  presented 
his  data. 

One  of  the  most  complete  and  thorough  experimental  investigations  of  the  turbulent 
boundary  layer  on  a  sharp  cone  at  incidence  under  supersonic  conditions  has  been  reported 
by  Rainbird  (Ref.  64).  The  experiments  were  performed  in  the  Canadian  National 
Aeronautical  Establishment  5-ft  intermittent  blowdown  wind  tunnel  at  Mach  numbers  of 
1.80  and  4.25  for  moderate  relative  incidence  (a/5v  «  1.2)  under  high  Reynolds  number 
conditions  with  a  relatively  hot  wall  (T* /T0 >so  =  1.0)  on  a  1 2. 5-dcg  semi-vertex-angle  sharp 
cone.  Complete  flow  conditions  are  given  in  Table  I  with  the  corresponding  windw;ard 
streamline  inviscid-edge  conditions  tabulated  in  Table  11.  Presented  in  Fig.  16  arc 
comparisons  of  calculated  local  Mach  number  profiles  across  the  turbulent  boundary  layer 
at  the  body  station  x/L  =  0.85  relative  to  experimental  measurements  (via  a  pitot  pressure 
probe)  taken  by  Rainbird  for  the  1.80  Mach  number  condition.  In  general,  the  agreement 
between  theory  and  experiment  is  reasonably  good  for  both  the  zero  angle  of  attack  and 
15.78-deg  angle-of-attack  conditions.  Note  that  the  experimentally  determined  turbulent 
boundary-layer  thickness  is  slightly  larger  than  predicted  by  the  present  theory.  Also  shown 
in  Fig.  16  are  rough  indications  of  the  physical  boundaries  for  the  inner  and  outer  regions 
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of  the  turbulent  boundary  layer  as  discussed  by  Bradshaw  (Ref.  67)  in  his  recent  review. 
The  location  of  the  inner-outer  boundary  at  y/5  «  0.2  can  easily  be  derived  from  Eq. 
(30)  by  equating  the  inner  and  outer  mixing  lengths  to  yield 


J-  =  — 

n  k* 


0.09 

0.435 


0.206 


OOl) 


where  yg  is  the  value  of  y  at  the  point  where  the  velocity  in  the  boundary  layer  is  equal 
to  0.99  of  the  velocity  at  the  boundary-layer  outer  edge,  and  hence  yg  «  8  where  6 
is  the  boundary-layer  thickness  defined  as  the  normal  distance  from  the  surface  where 
the  boundary-layer  velocity  equals  0.995  of  the  outer-edge  velocity. 


In  order  to  determine  the  local  Mach  number  distribution,  Rainbird  measured  the 
pitot  pressure  across  the  boundary  layer  which,  when  coupled  with  the  assumption  that 
the  static  pressure  remains  constant  through  the  boundary  layer  at  a  value  corresponding 
to  the  surface  static  pressure,  yields  the  local  Mach  number  distribution.  Rainbird  then 
assumed  the  total  temperature  to  remain  constant  across  the  boundary  layer  which  allows 
the  distribution  of  mean  velocity  to  be  calculated  through  the  boundary  layer.  With  this 
information,  integrated  boundary-layer  parameters,  such  as  the  displacement  thickness, 
which  depend  on  the  actual  boundary-layer  structure  can  be  determined  based  on  the 
experimental  results.  Figure  17  shows  a  comparison  of  Rainbird's  experimentally 
determined  (in  the  above-discussed  manner)  boundary-layer  parameters  at  the  body  station 
x/L  =  0.85  relative  to  calculated  distributions  from  the  present  windward  ray 
boundary-layer  analysis  for  both  Mach  number  conditions.  Note  that  the  present  theory 
underpredicts  the  boundary-layer  thickness  for  both  Mach  numbers  at  both  zero  and 
nonzero  angles  of  attack.  However,  the  calculated  two-dimensional  boundary-layer 
displacement  thicknesses  are  in  excellent  agreement  with  the  experimentally  determined 
values  by  Rainbird  for  all  cases.  It  should  be  noted  here  that  Rainbird  used  the 
two-dimensional  form  of  the  displacement  thickness  integral  to  evaluate  his  experimentally 
determined  value.  As  can  be  seen  from  Fig.  17,  the  three-dimensional  displacement 
thickness  is  slightly  smaller  than  the  two-dimensional  value  which  is  similar  to  the  results 
presented  previously  in  Fig  14  under  hypersonic  conditions. 


A  recent  note  by  Copper  and  Shaw  (Ref.  68)  was  concerned  with  presentation  of 
a  correlation  for  windward  ray  turbulent  boundary-layer  thickness  as  a  function  of  cone 
half-angle  and  angle  of  attack.  Based  on  examination  of.  a  large  amount  of  experimental 
data.  Copper  and  Shaw  recommend  the  simple  expression 


t?n[g(a^0)]  =  -2  +  2  exp[—  0.3a/8v]  (102) 

as  an  adequate  fit  to  the  data  for  engineering  purposes.  The  reader  is  referred  to  their 
Note  for  details  of  how  the  above  expression  is  to  be  modified  for  surface  mass  injection 
effects,  if  present.  Presented  in  Fig.  18  are  selected  results  comparing  the  present  analytical 
analysis  for  fully  turbulent  boundary  layers  with  the  Copper  and  Shaw  correlation.  In 
general,  good  agreement  is  observed.  For  the  case  of  a  laminar  boundary  layer.  Fig.  18 
shows  that  the  boundary-layer  thickness  will  decrease  more  rapidly  than  for  a  turbulent 
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boundary  layer  as  the  angle  of  attack  increases,  which  is  related  to  the  more  severe  crossflow 
influence  on  laminar  boundary  layers  relative  to  turbulent  boundary  layers  as  discussed 
previously. 


SECTION  IV 

CONCLUDING  SUMMARY 

The  present  investigation  was  devoted  to  the  formulation  and  application  of  a 
windward  plane  of  symmetry  laminar,  transitional,  and  turbulent  boundary-layer  analysis 
for  sharp  cones  at  incidence  in  a  supersonic  or  hypersonic  flow.  The  governing 
boundary-layer  equations  in  the  plane  of  symmetry  were  numerically  integrated  on  a  digital 
computer  using  an  implicit  finite-difference  technique  which  marched  along  the  windward 
streamline  starting  at  the  apex  of  the  cone  with  a  laminar  similar  solution.  A  so-called 
invariant  model  of  turbulence  was  used  in  a  two-layer  eddy  viscosity-mixing  length  approach 
for  calculation  of  the  turbulent  boundary  layer  in  conjunction  with  an  intermittency  factor 
treatment  of  the  transition  zone.  Inviscid  edge  conditions  were  input  to  the  boundary-layer 
analysis  from  a  previously  documented  inviscid  analysis  of  the  conical  flow  field  about 
a  sharp  cone  at  incidence  in  a  supersonic  or  hypersonic  flow. 

In  general,  the  agreement  between  the  present  theory  and  experimental  data  for 
measured  surface  heat  transfer  along  the  windward  streamline  was  good  to  excellent  for 
both  laminar,  transitional,  and  turbulent  boundary  layers  under  a  wide  range  of  flow  and 
angle-of-attack  conditions.  Comparison  of  calculated  and  measured  local  Mach  number 
profiles  in  a  fully  turbulent  boundary  layer  revealed  reasonable  agreement  with  the 
turbulent  boundary-layer  thickness  slightly  underpredicted  by  the  present  theory  relative 
to  experimental  measurements.  Integrated  boundary-layer  displacement  thickness 
calculations  were  found  to  compare  quite  well  with  experimentally  determined  values  using 
a  common  definition  of  the  displacement  thickness. 

The  effects  of  three-dimensional  crossflow  (outflow)  along  the  windward  plane  of 
symmetry  were  shown  to  be  significantly  stronger  in  a  laminar  boundary  layer  relative 
to  a  turbulent  boundary  layer.  This  finding  was  used  in  conjunction  with  an  "effective 
cone"  concept  to  show  that  windward  ray  heating  distributions  for  a  turbulent  boundary 
layer  could  be  computed  satisfactorily  (for  the  geometry  and  flow  conditions  under 
examination)  using  an  "effective  cone"  of  semivertex  angle  equal  to  the  physical  body 
semivertex  angle  plus  the  physical  angle  of  attack  in  a  two-dimensional  zero  angle-of-attack 
analysis.  Application  of  this  effective  cone  concept  to  the  windward  streamline  laminar 
boundary  layer  resulted  in  a  20-  to  30-percent  underprediction  of  heat  transfer  relative 
to  the  windward  plane  of  symmetry  analysis.  In  general,  smaller  crossflow  (outflow)  effects 
on  the  windward  streamline  boundary  layer  can  be  expected  for  turbulent  layers  as 
compared  with  laminar  layers  subject  to  the  same  boundary  conditions.  It  is  important 
to  realize  and  acknowledge  that  "equivalent  cone"  techniques  do  not  yield  uniformly  valid 
approximations  but  may  provide  acceptable  results  in  some  cases  (such  as  the  present 
plane  of  symmetry  turbulent  boundary  layer).  Caution  is  highly  recommended  in  the  use 
of  "equivalent  cone"  techniques  under  flow  and  angle-of-attack  conditions  where  their 
applicability  has  not  been  properly  assessed. 
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5.0-deg  Half-Angle  Sharp  Cone 
Mqj  =  12.0,  Re^/ft  -  9.5  x  105 
Vo,  oo  =  °- 14-  L  =  0.50ft 


Fig.  2  Comparison  of  Present  Results  with  Experimental  Data  from  Wittliff 
and  Wilson  (Ref.  58) 
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10.0-deg  Half-Angle  Sharp  Cone 
Mqj-8.0,  Re^/ft  =  1.3  x  106 

VTo,co  =  0-40'  L  =  0.333  ft 
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15.0-deg  Half-Angle  Sharp  Cone 
Mj,- 10.4,  Remm  -  3. 7  x  llj5 
00fl5r  TA®-°  °’  L- 0.271  ft,  x/L  -1.0 
j  Laminar  Boundary  Layer 
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10.0-deg  Half-Angle  Sharp  Cone 
Mq,  -  6.86,  Re^/ft  -4.6x  106 

V^co"054-  L"  1015  ft 
Laminar  Boundary  Layer 
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Experimental  Data  from 
Fig.  12b  of  Ref.  61 


Present  Theory 
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■  8.0 


4 


2 

1-3 

8 

6 

4 


2 


6.0-deg  Half -Angle  Sharp  Cone 
Mg,  *  10. 17,  Re^/ft  - 1.6  x  100 

Tv/ro, 00 =  °- 25-  L- 3.99ft 
Natural  Transition 


Present  Theory 


I - 1 - ■  i _ i _ I _ i _ i _ i  i _ l 

0  0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9  1.0 

x/L 


Fig.  6  Comparison  of  Present  Results  for  Natural  Transition  with  Experimental  Data  from 
McCauley,  Saydah,  and  Bueche  (Ref.  62) 
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6.0-deg  Half-Angle  Sharp  Cone 
Mg,  -  10. 17,  Re^/ft  - 1.6  x  1(£ 
yro<D-0.25,  L  =  3.99  ft 
Roughness- Induced  Transition 
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a.  Streamwise  Velocity  Ratio 
Fig.  8  Calculated  Laminar  Boundary- Layer  Profiles 
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y/L 


0 


Flow  Conditions  of  Fig.  6 
Body  Location  x/L  -  0.40 
Laminar  Boundary  Layer 
L- 3.99  ft 


c.  Static  Temperature  Ratio 
Fig.  8  Continued 
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d.  Total  Temperature  Ratio 
Fig.  8  Continued 
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e.  Local  Mach  Number 
Fig.  8  Continued 
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f.  Pitot  Pressure  Ratio 
Fig.  8  Concluded 
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a.  Streamwise  Velocity  Ratio 
Fig.  9  Calculated  Turbulent  Boundary- Layer  Profiles 
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b.  Crossflow  Velocity  Ratio 

Fig.  9  Continued 
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c.  Static  Temperature  Ratio 
Fig.  9  Continued 
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d.  Total  Temperature  Ratio 
Fig.  9  Continued 
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e.  Local  Mach  Number 
Fig.  9  Continued 
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f.  Pitot  Pressure  Ratio 
Fig.  9  Concluded 
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b.  Crossflow  Velocity  Ratio 
Fig.  10  Continued 


61 


AEDC-TR-71-235 


c.  Static  Temperature  Ratio 
Fig.  10  Concluded 
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a.  Local  Heat  Transfer 

Fig.  11  Comparison  of  Present  Results  with  "Effective  Cone"  Calculations 
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b.  Local  Skin  Friction 
Fig.  11  Concluded 
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6v  +  a,  deg 


6V  -  6. 0  deg,  a  4  0  deg 
Present  Theory 
6V  *  8. 0  deg,  10. 0  deg 
a  •  Odeg 

"Effective  Cone"  Theory 


a.  Local  Heat  Transfer 

Fig.  12  Summary  of  "Effective  Cone"  Results  Based  on  Fig.  11 


65 


AEDC-TR-71-235 


6V  +  a,  deg 


b.  Local  Skin  Friction 
Fig.  12  Concluded 
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6V  -6.0 deg,  (HO  deg 
Present  Theory 
6V  ■  8.0  deg,  10.0  deg 
a  a Odeg 

"Effective  Cone"  Theory 
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6.0-deg  Half -Angle  Sharp  Cone 


Fig.  13  Calculated  Boundary- Layer  Thickness  and  Displacement  Thickness 
Distributions 
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Fig.  14  Comparison  of  Two-  and  Three-Dimensional  Displacement  Thickness 
Distributions 
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a.  a  =  0  deg 

Fig.  16  Comparison  of  Calculated  Local  Mach  Number  Profiles  with  Experimental 
Measurements  from  Rainbird  (Ref.  64) 
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b.  a  =  15.78  deg 
Fig.  16  Concluded 
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12.5-deg  Half -Angle  Sharp  Cone 
Mq,  =  1.80,  ReQj/ft  =  7.4  x  10^ 

Woo'10'  l- =  3.46 ft 
Turbulent  Boundary  Layer 

Mf'=0.995/L  ) 

- 6J/L  /  Present 
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b.  M„  =  1.80 
Fig.  17  Concluded 
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-Curve  Fit  to  Turbulent  Experimental  Data 
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Fig.  18  Comparison  of  Calculated  Boundary- Layer  Thicknesses  with  Correlation  to  Experimental 
Data  as  Proposed  by  Copper  and  Shad  (Ref.  68) 
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TABLE  I 

FLOW  CONDITIONS 


Ref. 

6v,deg 

Rejft 

Tw/To,- 

Po  ,»)Psia 

To,-,°R 

L,in. 

58 

5.0 

12.0 

9.5  x 

10s 

0.14 

4200.0 

3600.0 

6.0 

59 

10.0 

8.0 

1.3  x 

106 

0.40 

259.0 

1360.0 

4.0 

60  ‘ 

15.0 

10.4 

3.7  x 

106 

0.23 

3830.0 

2260.0 

3.25 

61’ 

10.0 

6.86 

4.6  x 

106 

0.54 

452.0 

1020.0 

12.2 

62 

6.0 

10.17 

1.6  x 

106 

0.25 

1310.0 

1875.0 

47.8 

63 

8.0 

10.2 

2.1  x 

10« 

0.27 

1710.0 

1890.0 

35.9 

64 

12.50 

4.25 

1.5  x 

107 

1.0 

166.0 

530.0 

41.6 

64 

12.50 

1.80 

7.4  x 

106 

1.0 

25.0 

530.0 

41.6 
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TABLE  II 

WINDWARD  STREAMLINE  INVISCID  EDGE  CONDITIONS 


Ref. 

5v,deg 

a,  deg 

P/P- 

Te/T. 

Me 

-x 

58 

5.0 

0.0 

12.0 

2.817  . 

1.371 

10.180' 

0.0  ' 

58 

5.0 

1.0 

12.0 

3.467 

1.481 

9.777 

1.556  x  10-3' 

58 

5.0 

2.0 

12.0 

4.232 

1.608 

9.364 

3.283  x  10-3 

58 

5.0 

3.0 

12.0 

5.114 

1.751 

8.950 

5.124  x  10-3 

58 

5.0 

4.0 

12.0 

6.109 

1.911 

8.543 

7.080  x  10-3 

58 

5.0 

5.0 

12.0 

7.220 

2.089 

8.145 

9.191  x  10-3 

59 

10.0 

0.0 

8.0 

4.068 

1.584 

6.210 

0.0 

59 

10.0 

2.0 

8.0 

5.216 

1.771 

5.827 

5.143  X  10*3; 

59 

10.0 

4.0 

8.0 

6.557 

1.987 

5.452 

1.104  x  lO-2 

59 

10.0 

6.0 

8.0 

8.084 

2.231 

5.092 

1.768  x  lO*2 

59 

10.0 

8.0 

8.0 

9.789' 

2.503 

4.751 

2.521  x  lO2 

59 

10.0 

10.0 

8.0 

11.668 

2.801 

4.431 

3.384  x  10  2 

60 

15.0 

0.0 

10.4 

11.886 

2.858 

5.882 

0.0  * 

60 

15.0 

3.0 

10.4 

16.203 

3.546 

5.187 

7.935  x  lO*3 

60 

15.0 

6.0 

10.4 

21.164 

4.335 

4.594 

1.829  x  lO-2 

60 

15.0 

9.0 

10.4 

26.706 

5.214 

4.087 

3.145  x  lO*2 

60 

15.0 

12.0 

10.4 

32.759 

6.172 

3.652 

4.790  x  10-2 

61 

10.0 

0.0 

6.86 

3.309 

1.456 

5.546 

0.0 

61 

10.0 

2.0 

6.86 

4.147 

1.595 

5.258 

5.803  x  10-3 

61 

10.0 

4.0 

6.86 

5.125 

1.754 

4.968 

1.241  x  lO2 

61 

10.0 

6.0 

6.86 

6.241 

1.934 

4.682 

1.973  x  lO-2 

.  61 

10.0 

8.0 

6.86 

7.488 

2.133 

4.405 

2.789  x  lO  2 

62 

6.0 

0.0 

10.17 

2.872 

1.381 

8.575 

0.0 

62 

6.0 

2.0 

10.17 

4.056 

1.579 

7.980 

3.934  x  10-3 

62 

6.0 

4.0 

10.17 

5.567 

1.824 

7.378 

8.512  x  lO  3 

63 

8.0 

4.0 

10.20 

7.599 

2.155 

6.753 

9.003  x  lO*3 

64 

12.50 

0.0 

4.25 

2.454 

1.307 

3.556 

0.0 

64 

12.50 

15.625 

4.25 

6.686 

2.000 

2.556 

1.036  x  10-1 

64 

12.50 

0.0 

1.80 

1.360 

1.092 

1.596 

0.0 

64 

12.50 

15.78 

1.80 

2.066 

1.234 

1.295 

2.558  x  10-1 
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APPENDIX  III 

IMPLICIT  FINITE-DIFFERENCE 
SOLUTION  OF  GOVERNING 
BOUNDARY-LAYER  EQUATIONS 

For  the  special  case  of  a  sharp  cone  at  zero  to  moderate  angle  of  attack  (0  =  a 
«  1.2  6V)  in  a  supersonic  or  hypersonic  flow,  the  inviscid  flow  field  is  conical  in  character 
so.  that  all  the  inviscid  flow  parameters  (p,  Te,  pe,  Ue,  Me,  3We/30,  32p/3 02)  remain 
constant  along  rays  from  the  apex  of  the  cone.  Under  this  restriction  the  relationship 
between  £  and  x  is  found  by  integration  of  Eq.  (42)  with  r  =  x  sin  5V  from  geometry 
which  results  in 

t  =  sin2SvX3  (HI-1) 

This  definition  of  £  leads  to  a  simple  expression  for  the  G  function  defined  by  Eq.  (62), 
namely 

r  _  _ 2 _ 

U  “  3  sin  Sv  (III-2) 

i  y 


Noting  that  dr/dx  =  sin  5V  geometry  and  (3=0  due  to  the  inviscid  conical  flow,  the 
governing  (Jane  of  symmetry  boundary-layer  equations  (51),  (52),  and  (53)  in  transformed 
coordinates  reduce  to  the  following  simpler  set: 

x-MOMENTUM 


*  Ef  -  1  +  Ge]'"  * 

^-MOMENTUM 

£V'-  +  t  +  Gg]g"-  G[(g')2  +  f'g'  sin  +xe]  = 

ENERGY 

e "  +  (Q  +  f  +  G«]0'  +  (y-i>Me*ttf">2  = 


(III-3) 


(III-4) 


(III-5) 


The  boundary  conditions  remain  unchanged  from  Eqs.  (65)  through  (68)  which  are  repeated 
below  for  completeness. 


MOMENTUM 

f(£,  T)  =  0)  =  0 
f'<£  T7  =  0)  =  0 


(III-6) 

(III-7) 
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l'(f, ,) .  1 

rj^oc 

(III-8) 

g(^,  ?7  =  0)  =  0 

(IH-9) 

g'(£  17  =  0)  =  0 

lim  ,  „  . 

g'(£  =  g'e 

7^-»oo  e 

(III- 10) 

(III-ll) 

ENERGY 

0(£  77  =  0)  =  8W  [0 

(III- 12) 

lim 

Troa  «£  v)  =  i 

(III-13) 

Following  the  approach  by  Blottner  (Ref.  42)  and  Davis  (Ref.  69)  the  momentum 
and  energy  equations  (HI-3,  -4,  -5)  are  rewritten  in  "standard"  form  for  a  parabolic  partial 
differential  equation  as 

<?2w  <?w  aw 

+  “1  +  Q2W  +  a3  +  a4^  = 

:  0  (III- 14) 

where  W  =  f*  for  the  x-momentum  equation,  W  =  g'  for  the  0-momentum  equation,  and 
W  =  6  for  the  energy  equation.  Using  Eqs.  (III-3,  -4,  -5),  one  finds  the  coefficients  a\ 
through  04  in  linearized  form: 

x-MOMENTUM 

dl*  .  r 
—  -  i  +  Qg  +  2^ 

ai  e* 

(III- 15) 

a2  =  0 

(III- 16) 

o 

II 

CO 

(III-17) 

-2ff' 

“4  £* 

(III- 18) 

^-MOMENTUM 

dl"  ,  r  dl 

dj  +  f  +  Gg+  2^ 

*'  = 

(III-19) 

2f ' 

-Ce'-T- 
“2  -  .* 

(III-20) 
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-GX0 

“3  =  * 

3 

(III-21) 

-2£f' 

“4  -  g* 

(III-22) 

ENERGY 

(t) 

(III-23) 

“2  =  0 

(III-24) 

(y-l)MeZ  £*(f")2 

“3_  (t) 

(111-25) 

-2ff' 

‘  f) 

(III-26) 

The  rj  derivatives  in  Eq.  (III-14)  are  replaced  with  finite-difference  quotients  which 
allow  variable  grid  spacing  in  the  r?  direction  so  as  to  concentrate  grid  points  in  the  region 
near  the  body  surface  where  the  dependent  variables  change  most  rapidly.  The  derivative 
in  the  £  direction  in  Eq.  (Ill- 14)  is  handled  in  the  usual  manner  as  a  two-point  backward 
difference  between  points  (m+l,n)  and  (m,n),  whereas  all  rj  derivatives  are  evaluated  at 

point  (m+l,n)  according  to  the  grid  mesh  shown  below. 


Station 
ra  m+1 
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The  solution  is  assumed  to  be  known  at  point  (m,n)  and  unknown  at  point  (m+l,n) 
so  that  the  finite-difference  scheme  to  be  constructed  will  be  implicit  in  nature.  The 
finite-difference  replacements  for  the  derivatives  are  as  follows  (see  Appendix  IV  for 
derivation): 


where 


[d2\ 

?1 

2[W  .  ,  -r  KW  -  ( 1 4-  K)  W  ]  .  , 

n  +■  1  n—  1  n  m  +  1 

(III-27) 

1 

D„ 

L  v 

*  m  +  X,  n 

2 

few 

l*.+  i  - 

(III-28) 

m  +  1,  n 

<?w 

w  .  ,  -  w 

m  +  1,  n  m ,  n 

(III-29) 

5F 

m+  1,  n 

Di  = 

(7n+l  "  7n)  +  K^n-^n-l) 

(III-30) 

°2  = 

(7n+  1  “  7n>2  +  K(»?n  -  7n-l)2 

(III-31) 

K  = 

7n  +  1  -  7n  .  . 

- (constant; 

(III-32) 

7„-7n-i 

n  = 

tj2  -  Jj]  (constant) 

(III-33) 

The  finite-difference  form  of  Eq.  (Ill- 14)  becomes,  upon  substitution  of  Eqs.  (III-27,  -28, 
-29), 


m+  1,  n+  1 


Bn  Wm+l,n 


Cn  Wm+1,  n—  1 


=  R. 


(III-34) 


where 


.  2  al 

n  =  D2  +  Dj 


B_  = 


-201 +  K)  a  j(l  —  K 2)  a4 

- -  +  “2  +  Af 


C  =  ^  - 


2K  K^a , 


°2  D1 


R„  =  -  a,  + 


a  a  w 

4  m ,  n 


AT 


(III-35) 


(III-36) 

(III-37) 

(III-38) 


In  order  . for  Eq.  (III-34)  to  be  linear,  the  coefficients  An,  Bn,  Cn,  and  Rn  must  be  treated 
as  known  quantities  at  point  n;  more  will  follow  on  this  subject  later.  The  important 
point  is  that  Eq.  (III-34)  represents  a  set  of  simultaneous  linear  algebraic  equations  under 
this  restriction. 
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Since  the  simultaneous  linear  algebraic  equations  resulting  from  Eq.  (III-34)  are  df 
a  special  form  (tridiagonal),  an  efficient  method  of  solution  on  a  digital  computer  is 
available  from  Richtmyer  and  Morton  (Ref.  70,  pp.  198-201,  274-282).  For  this  procedute 
the  boundary  conditions  at  the  wall  (n=l)  and  the  outer  edge  (n=N)  must  have  specified 
values  Wm  +  1>1  and  Wm  +  i(N.  Due  to  the  special  form  of  the  Eq.  (III-34),  the  relation 

Wm+1,  n  =  En  *«,+  !,  n+1  +  en’  2<n^N-l  (IH-39) 

exists  where 

E„  =  ^  (III-40) 

B2 


e2 


r2  -  c2  w 


m+  1,  1 


(III-41) 


-A. 


E„  " 


Bn  +  Cn  En  , 
n  n  n  —  i 


R- — C_  e  1 , 

n  n  n  —  i 

B  +  C  E  , 

n  n  n—  1 


3  <  n  <  N-l 


(1 11-42) 

(III-43) 


The  quantities  En  and  en  are  computed  from  Eqs.  (III-40)  through  (III-43)  starting  with 
n  =  2  and  progressing  to  n  =  N-l.  The  solution  Wm  +  i,n  is  then  obtained  by  evaluating 
Eq.  (III-39)  from  n  =  N-l  to  n=2.  Knowing  the  distribution  of  f  and  0  across  the  boundary 
layer  from  the  above  procedure,  the  transformed  stream  function  f  is  evaluated  from 


f  (£-?)  =  P  f  '(£>?)  dr, 

o 


(III-44) 


where  the  integral  is  numerically  integrated  using  the  well-known  trapezoidal  rule,  viz., 


-iE  7-^)  + 

i=  2 


(III-45) 


with 


Di  = 


7i-T?i_  i 


(III-46) 


The  transformed  stream  function,  g,  is  evaluated  in  a  similar  manner.  Inversion  from  the 
transformed  (£,tj)  plane  to  the  physical  (x,y)  plane  is,  from  Eq.  (85), 


y  = 


#  ru,  *, 

PeUer  J  p 


(III-47) 
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where 


(III-48) 


due  to  the  constancy  of  static  pressure  across  the  boundary  layer.  Again  using  the 
trapezoidal  rule  method  of  numerical  integration  yields 

prr  n 

yn  =  TTTtE  +  ^7  =  7i_i)lDi  (III-49) 

*  e  e 

i=  2 


with  Dj  given  by  Eq.  (III-46).  The  relationship  between  \  and  x  for  sharp  cone  flows 
is  given  by  Eq.  (III-l)  which  results  in 


(III-50) 


For  the  present  sharp  cone  study,  Eq.  (III-50)  was  used  in  conjunction  with  an  x-step 
size  division  Ax  =  0.01  L,  i.e.,  100  stations  spaced  equally  along  the  physical  body. 


The  mathematical  basis  of  the  above  tridiagonal  matrix  procedure  applied  to  the 
solution  of  boundary-layer  problems  is  due  to  Fliigge-Lotz  and  Blottner  (Ref.  71).  The 
present  application  differs  from  their  original  work  in  one  important  aspect,  the  linearized 
difference  equations  herein  are  uncoupled  and  solved  separately.  In  Fliigge-Lotz  and 
Blottner 's  approach,  the  difference  equations  remain  coupled  and  require  additional 
machine  storage  and  manipulations  for  solution.  With  the  present  uncoupled  approach, 
the  difference  equations  are  iterated  to  convergence  at  each  station  along  the  body  so 
that  one  must  pay  the  price  of  iteration:  this  procedure  is  how  the  linearizing  coefficients 
a.i  through  04  are  evaluated  at  each  station  using  the  results  of  the  previous  iteration. 
Iteration  at  a  given  station  is  continued  until  successive  values  of  f',  g',  and  d  differ  by 
less  than  0.10  percent  from  the  corresponding  value  of  the  preceding  cycle.  Typically 
about  three  iterations  per  station  are  required  for  fully  turbulent  flow;  only  one  to  two 
iterations  per  station  are  necessary  for  laminar  flow.  By  use  of  the  above  procedure  the 
final  solution  obtained  at  each  station  is  exact  in  the  sense  that  it  represents  a  converged 
iterated  solution  to  the  governing  nonlinear  partial  differential  equations  written  in 
finite-difference  form. 


The  variable  grid  mesh  used  in  the  present  work  is  taken  from  Smith  and  Cebeci 
(Ref.  2).  The  various  constants  used  herein  are  as  follows: 

N  =  85 
K  =  1.063 

a  =  0.010 

Experience  with  varying  these  constants  and  observing  their  influence  on  the  resultant 
numerical  solution  has  indicated  that  the  above  choices  are  adequate  under  the  present 
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flow  conditions.  These  values  may  not  be  satisfactory,  however,  for  other  body  geometries 
and  flow  conditions  so  that  the  influence  of  the  variable  grid  mesh  constants  should  be 
ascertained  for  each  new  investigation.  Provisions  are  included  in  the  present  program  for 
addition  of  points  (if  necessary)  to  the  solution  as  the  body  is  traversed.  This  may  be 
required  in  order  to  prevent  numerical  supression  of  boundary-layer  growth,  especially 
in  the  case  of  transitional  and  turbulent  boundary  layers. 

Evaluation  of  heat-transfer  rate  and  shear  stress  at  the  cone  surface  requires  numerical 
determination  of  the  wall  derivatives  0'(£,tj=O)  and  f"(£,rpO)  since 

4*  «  0  '(£  7  =  0) 

.  •  rw  «  f"(£  7  =  0) 

from  Eqs.  (82)  and  (84).  Due  to  the  variable  grid  spacing  in  the  17-direction  discussed 
previously,  the  wall  derivatives  are  evaluated  in  the  present  work  by  application  of  the 
three-point  Lagrangian  interpolation  formula  (Ref.  72,  pp.  71-77)  evaluated  according  to 
the  following  diagram 


which  yields,  at  station  m+1, 

f"(7  =  0)  =  n[B1(2+K)f',i  +  +  B3Vv]  (HI-51) 

0'(7  =  O)  =  il[B.(2  +  K)0  -  B2(1+K)0_  -  B30_  ]  (III-52) 

1  Vi  z  72  *  v3 
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where 


-l 

q2(uk) 


(III-53i) 


(III-54) 


and 


-1 

(III-55) 

Q2KU  +  K) 

=  f'(7?  =  i?5) 

(III-56) 

=  6(r}=T)J 

(III-57) 

with  i  =  1,  2,  3.  The  parameters  K  and  S2  are  defined  by  Eqs.  (III-32)  and  (III-33), 
respectively.  The  use  of  the  three-point  Lagrangian  interpolation  formula  for  the  wall 
derivatives  is  consistent  with  the  formulation  of  the  variable  grid  spacing  derivatives 
discussed  in  Appendix  IV  and  follows  the  work  of  Smith  and  Cebeci  (Ref.  2,  Appendix 
C). 


At  the  apex  of  a  sharp  cone  where  £  =  0,  the  governing  boundary-layer  equations 
(III-3),  (III-4),  and  (III-5)  reduce  to  the  following  ordinary,  nonlinear,  differential 
equations: 

x-MOMENTUM 

*V"  +  +  f  +  Gg]f"  =  0  (III-58) 

0-MOMENTUM 

£* g'"  +  [|^  +  f  +  GgJg"  -  G[(g')2  +  f'g'sin5v  +  X6]  =  0  (III-59) 


ENERGY 


(pt)0''  +  [^(0  +  f  '  Gg]0'  +  (y-1)  Me£^r')2  =  0  (HI-60) 

with  the  boundary  conditions  (III-6)  through  (III-13).  In  order  to  obtain  starting  profiles 
for  the  finite-difference  scheme  to  march  downstream,  Eqs.  (III-58),  (III-59),  and  (III-60) 
are  solved  by  using  the  tridiagonal  matrix  procedure  described  previously  in  the  following 
manner.  At  the  apex  of  the  sharp  cone,  £  =  0  so  that  <24  =  0  in  both  of  the  momentum 
equations  as  well  as  the  energy  equation;  in  addition,  the  term  containing  £  in  a,\  of 
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both  the  momentum  and  energy  equations  vanishes.  Initial  guesses  for  f,  g\  and  Q  are 
input  to  the  program  as 


f'j  =  1  -  exp(-7/) 

(III-61) 

«I  =  6'ef'l 

(III-62) 

*1  =  K  +  b-Wl 

(III-63) 

where  the  subscript  I  denotes  the  initial  approximation.  The  equations  are  then  iterated 
to  convergence  in  the  same  manner  described  previously;  an  averaging  scheme  is  used  to 
speed  convergence.  Typically  about  20  to  25  iterations  are  required  to  generate  a  converged 
initial  solution. 

By  use  of  the  above  procedures,  the  numerical  solution  of  any  two-point 
boundary-value  problem  governed  by  either  linear  or  nonlinear  ordinary  differential 
equations  as  well  as  sets  of  coupled  parabolic  partial  differential  equations,  either  linear 
or  nonlinear,  is  reduced  to  subroutine  status  on  a  digital  computer  in  that  only  the 
coefficients  a\  through  04  must  be  defined  in  conjunction  with  the  required  boundary 
conditions  for  each  new  problem.  Much  application  of  these  procedures  has  been  found 
in  the  von  Kirmdn  Facility  of  the  Arnold  Engineering  Development  Center  on  many 
different  problems  involving  both  nonlinear  ordinary  differential  equations  as  well  as 
parabolic  partial  differential  equations.  The  thin  viscous  shock  layer  analyses  by  Adams 
(Refs.  73  and  74)  are  good  examples.  In  these  works  eight  simultaneous,  nonlinear,  ordinary 
differentia]  equations  governing  momentum,  energy,  and  species  conservation  were  solved 
using  the  iterative  tridiagonal  matrix  method  presented  above.  The  problem  was  made 
even  more  difficult  in  that  chemical  nonequilibrium  effects  were  included  in  the  analysis. 
Based  on  experience  with  analyses  of  this  type,  the  author  highly  recommends  use  of 
the  iterative  tridiagonal  matrix  approach  where  applicable. 


85 


AEDC-TR-71-235 


APPENDIX  IV 

VARIABLE  GRID  DIFFERENTIATION  FORMULAS 

The  first  and  second  derivative  formulas  for  the  variable  grid  system  discussed  in 
Appendix  III  can  be  most  easily  derived  from  consideration  of  the  three-point  Lagrangian 
interpolation  polynomial  (Ref.  72,  pp.  71-77)  evaluated  according  to  the  following  diagram:' 


which  yields,  for  the  function  W(tj)  in  the  range  r/a.i  <  tj  <  Tja+i, 

n+  1 

W(,)  =V  Li(,)  w.  (7,j)  *  E(,)  (IV-1) 

i  =  n-l 


where 

n  + 1 
j=n  — 1 

W*  -  (,V-2) 

j  =  n  —  1 

j^k 


and  E(tj)  represents  the  error  of  the  interpolating  polynomial  as  discussed  in  Section  3.2 
of  Ref.  72.  Differentiating  Eq.  (IV-1)  yields 


i  =  n-l 


(IV-3) 
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d2W(7?) 

dr,2 


■s 


d^L;(7) 


w 


(IV-4) 


where  [dLi(q)]/dr)  and  [d2LiOj)]/dr?2  are  to  be  evaluated  from  differentiation  of  Eq.  (IV-2). 

The  present  application  is  concerned  with  evaluation  of  the  above  derivatives  at  the 
location  77  =  r?n.  Under  this  restriction,  Eq.  (IV-3)  can  be  expanded  to  read 


dW(r?n) 

dr? 


(,?n  “  Vn+l> 


(,?n-l-,?n)(,?n-l-,/n-l) 


~  +  {}?„  -  *?„+!>  „„  , 


(,?n-,?n-l) 


(j?n+]  -  I7n_])  (r?n+]  -  r?n) 


W(r?n.,) 


dE(r?n) 

d»? 


Introducing  the  nomenclature 


Ar?+  =  r?n+1  -  r?n 
and 

A»?-  =  7„  "  7n— 1 


(IV-5) 


(IV-6) 


(IV-7) 


Equation  (IV-5)  can  be  written  as,  by  forming  a  common  denominator  for  the  first  three 
terms  on  the  right-hand  side, 


dW(r?n) 

dr? 


-(A7?+)2W(r?n_i)  -  [(Ar?_)2  -  (At?+)2]  Wh?n)  +  (Ar?J2W(r?nH  x) 
[A77JA77J2]  +  (At?_(At?+)  ^ 


dE(r?  ) 

-  sr  (IV'8) 


or,  upon  division  of  both  numerator  and  denominator  of  the  first  term  on  the  right-hand 
side  by  (At?-)2, 


dW(r?n) 

dj? 


*  [  -  (a^)  JW(^  +  W(W 

7  77HY 

Ar?+  + 


+ 


dE(r?n) 

dr? 


(IV-9) 


Smith  and  Cebeci  (Ref.  2)  have  proposed  a  variable  grid  system  such  that  the  ratio 
of  lengths  of  any  two  adjacent  intervals  is  a  constant,  i.e., 


Ar?+  =  K  At? 


(IV-10) 
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where  the  value  of  K.  must  be  specified.  The  distance  from  the  body  surface  to  the  i-th 
grid  point  is  given  by  the  following  formula 

*■  -  q[ttt]  (IV-II) 

where  fl  is  the  length  of  the  first  step,  i.e., 

o  =  12  -  ’ll  (IV- 1 2) 

In  Eq.  (IV-12),  t?i  corresponds  to  the  grid  point  on  the  body  surface  (t?j  =  0)  and  172 
is  the  first  grid  point  away  from  the  wall.  Adopting  the  Smith  and  Cebeci  grid  system 
for  the  present  work  permits  Eq.  (IV-9)  to  be  written  as 


with 


dff(^)  -K2W(i?n_])  -  [1-K 2]W(»Jn)  +  W(7/n+1)  dE(7jn) 
~dr]  =  D~J  drj 


(IV- 13) 


D1  =  A,4  +  +  K2(i?n-l?n-l) 


(IV- 14) 


In  a  similar  manner,  Eq.  (IV-4)  can  be  expanded  to  read 


V 


2Ajj_W(7jn_1)  -  2(Ai/+ +  Ai7_)8{rjn)  +  2\r}_V(rfn+1)  d2E(?7n) 

A7+Aj?_(A7?+  +  A77J  +  drj2  CIV-1 5) 


which  becomes,  upon  division  of  the  numerator  and  denominator  of  the  first  term  on 
the  right-hand  side  by  Arj., 


V 


feHn-ii  -  (l  +  ^)w,n) +  w^,)| 

<H>2  -  (^J2 


Introducing  K  according  to  Eq.  (IV-10)  into  Eq.  (IV-16)  gives 


with 


(IV-16) 


d2W(77n)  2[KW(77n_j)  -  (1  +  K)ff(ijn)  +  W(jjn+1)]  d2E(7jn)  (IV-17) 

dTj2  °2  djj2 


D2  =  (Ar/+)2  +  UArjJ2  =  (77^1  -  7jn)2  +  Ui}n  -  r}n_})2  (IV-I8) 
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The  error  terms  can  be  shown  to  be  (see  Ref.  72,  pp.  108-113) 


dr] 


6 


d3W(n) 

—7T~  Ai4Ai?- 
dr) 


(IV- 19) 


d2E(?n) 

dr,2 


l 

12  dr,4 


(IV-20) 


where  tj,  is  some  7,-point  in  the  interval  r,n.i  ^  rj,  ^  r,n+i  so  that  theoretically,  by 
letting  A rj+  and  At,,  approach  zero,1  the  discretization  error  also  approaches  zero.  It  is 
implicitly  assumed  in  all  of  the  foregoing  analyses  that  the  function  W  and  its  derivatives 
are  single-valued,  finite,  and  continuous  functions  of  the  independent  variable  r,. 
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